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CHAPTER  I 


IHTRODUCTION  AND  SURVEY 

The  research  which  Is  the  subject  of  this  thesis  consists  of  (l)  the 
development  of  theories  and  kinetic  expressions  for  thermoluminescence 
from  macroscopic  crystals,  molecular  crystals  and  amorphous  solids, 

(2)  the  measurement  of  thermolxominescence  from  amino  acids,  polypeptides 
and  proteins  which  had  been  subjected  to  ultraviolet  irradiation,  (3)  the 
measurement  of  electron  paramagnetic  resonance  developed  in  these  materials 
as  a  result  of  ultraviolet  irradiation. 

Proteins,  polypeptides  and  their  basic  subunits,  amino  acids,  are 
ubiquitous  and  essential  constituents  of  all  biological  systems.  The 
polymers  serve  as  structural  units,  as  catalysts  of  biochemical  reactions, 
as  hormones  and  as  other  fxmctional  units.  Proteins  and  polypeptides  are 
composed  of  one  or  a  few  chains  of  mixed  polymers  in  which  the  various 
amino  acid  molecules  are  joined  together  with  peptide  bonds. ^  The  chains 
in  turn  are  held  together  with  covalent  or  hydrogen  bonds  and  possibly  by 
van  der  Waals  forces,^’® 

The  generalized  structural  formula  of  amino  acids  is  shown  below. 

R 

I 

HOOC  -  CH  - 

R  stands  for  residue  or  side  group.  There  are  about  twenty  different 
side  groups  found  in  naturally  occurring  amino  acids.  They  range  from 
R  =  H  in  glycine  to 

"  ■"  Co 

N 

1 

H 

^J.  L.  Oncley,  Rev.  Mod.  Phys.  30  (l959) • 

2a.  Rich,  Rev.  Mod.  Phys.  50  (1959). 
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in  tryptophan.  Polymerization  of  amino  acids  takes  place  with  the 


splitting  out  of 


HgO.  The  reaction  is 

HOOC  -  CH  -  MH.  +  HOOC  -  CH  - 

1  2  I  2 


H 

I 

HOOC  -  CH  -  N  -  C 


CH 

I 


-  NHg  +  HgO 


and  may  be  repeated  at  either  end  of  the  molecule.  and  Eg  may  be 

any  of  the  twenty  s Ide  groups . 

On  the  basis  of  x-ray  diffraction  and  other  studies  by  many  workers, 
it  is  now  believed  that  the  polypeptides  may  assume  three  different 
configurations.  These  are:  (l)  the  extended  zig-zag  or  pleated-sheet 
configuration  shown  in  Pig.  1,  (2)  the  Q!-hellx  configurations^  shown  in 
Pig.  2  and  (3)  the  random  coil.  The  pleated-sheet  configuration  is  found 
in  structural  or  filamentous  protein.  The  a-helix  and  the  random  coil 
configurations  are  both  found  in  many  globular  and  catalytic  proteins. 

Enzymes,  i.e.,  catalytic  proteins,  are  able  to  promote  a  wide  variety 
of  biochemical  reactions  without  being  used  up  in  the  process.  Of  parti¬ 
cular  interest  from  the  point  of  view  of  the  physics  Involved  are  the 
charge-transfer  reactions. 

Such  reactions  may  be  separated  into  two  groups  according  to  whether 
they  are  ordinary  collision  mediated  reactions  or  not.  A  collision 
mediated  reaction  is  one  in  which  two  molecid.es  moving  about  in  solution 
collide  and  form  an  excited  complex.  This  complex  then  breaks  apart  into 
molecules  different  from  the  initial  reactants,  A  common  example  is 

AHg  +  B  -  A  +  .BH'  +  H"^  , 


where  A  and  B  are  molecules  whose  exact  structures  are  of  no  concern 
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—(a)  A  single  polypeptide  chain,  (b)  Polypeptides  arranged  i 
pleated  sheet  configuration  [from  E.  E.  Marsh,  R.  B.  Corey, 
L.  Pauling,  Biochim.  et  Biophys,  Acta  1  (1955)]. 


In  charge  transfer  reactions  which  are  not  collision  mediated,  the 
charge  transfer  may  take  place  across  distances  of  from  10  to  50  Angstroms 
or  more.  Such  reactions  are  best  discussed  in  terms  of  examples,  even 
though  the  reaction  mechanisms  for  these  examples  are  not  completely 
understood. 

The  biological  electron  transport  system  consists  of  oxidation- 
reduction  enzymes  which  are  located  in  the  membrane  of  the  mitochondria. 
Electrons  move  throi;igh  this  sequence  of  enzymes  driving  the  formation  of 
energy-rich  compounds.  The  enzymes  pass  the  electron  to  one  another 
apparently  without  the  use  of  intermediate  shuttle  molecules.^  Further¬ 
more,  the  enzymes  are  restricted  to  certain  sites  in  or  on  the  membrane. 
They  may  or  may  not  have  complete  freedom  of  rotation  at  these  sites. 

Chance'^  has  proposed  two  mechanisms  for  charge  transfer  among  these 
enzymes.  In  one  mechanism,  freedom  of  rotation  of  the  molecules  is 
posttU-ated  to  permit  collision  mediated  charge  exchanges  between  the 
prosthetic  groups  of  neighboring  enzymes .  In  the  other  mechanism,  the 
molecules  are  assumed  to  be  fixed,  and  the  charges  move  through  them. 
These  two  mechanisms  are  schematically  indicated  in  Fig.  3-  In  these 
diagrams,  the  circles  represent  the  protein  part,  the  straight  line  the 
prosthetic  group  of  the  enzyme.  Arrows  in  the  upper  diagram  indicate  the 
notion  that  the  enzymes  rotate  until  the  prosthetic  groups  of  two  neigh¬ 
boring  enzymes  are  touching,  whereupon  charge  transfer  takes  place  from 
one  group  to  the  other.  In  the  lower  diagram,  the  arrows  indicate  a 
hypothetical  migration  of  charge  through  the  protein  from  one  prosthetic 
group  to  the  next.  These  diagrams  are  a  great  oversimplification  in  that 

^A.  L.  Lehninger,  Rev.  Mod.  Phys .  I36  (1959)- 

■^B.  Chance  and  G.  R.  Williams,  Advances  in  linzymol.  65  (l959)  • 


FIG.  3"“Schematic  representation  of  two  mechanisms 
of  charge  transfer  along  a  chain  of  fixed 
molecules  (adapted  from  Chance  and  Williams^). 
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they  leave  out  the  chemistry  known  to  be  associated  with  the  electron 
transfer.  Information  and  references  on  this  chemistry  may  be  found  in 
Lehninger's  article.^  Though  neither  of  these  mechanisms  has  been  sub¬ 
stantiated,  the  second  does  serve  as  an  example  of  a  charge  transfer  which 
takes  place  without  a  collision  of  two  molecules. 

A  second  example  is  also  the  work  of  Chance.®  He  has  studied  a  light- 
induced  intermolecular  oxidation-reduction  reaction  between  bacterial 
chlorophyll  and  a  cytochrome,  and  has  found  that  this  reaction  still 
proceeds  at  temperatures  as  low  as  80°  K.  Spectroscopic  data  show  that 
in  this  reaction,  the  cytochrome  is  oxidized  from  the  ferrous  to  the 
ferric  state.  The  quantum  requirement  for  this  reaction  is  two  photons 
per  transferred  charge.  The  photons  are  absorbed  by  chlorophyll  molecules. 
Since  only  three  percent  of  the  chlorophyl  molecules  are  contiguous  with 
cytochrome  molecules,  it  follows  that  the  energy  produced  when  any 
chlorophyll  molecule  absorbs  a  photon  can  effect  the  oxidation  over  a 
distance  large  compared  to  the  size  of  the  molecule.  And  since  the 
reaction  proceeds  at  80°  K,  a  bimolecular  collision  is  clearly  not  involved. 
Several  explanations  are  possible  as  to  how  the  oxidation  takes  place. 

(1)  The  photon  excites  an  electron  In  the  chlorophyll  and  it  or  its 
vacancy  migrates  to  the  site  of  a  cytochrome  where  oxidation  takes  place. 

(2)  The  excitation  caused  by  the  photon  migrates  as  an  exciton  to  the 
cytochrome  and  there  promotes  the  oxidation.  (3)  The  photon  may  initiate 
a  proton  migration®  which  culminates  in  the  oxidation.  An  Inductive 
resonance  between  an  excited  chlorophyll  and  a  cytochrome  is  unlikely 
because  the  emission  spectrum  of  chlorophyll  and  the  absorption  spectrum 

®B.  Chance  and  M.  Nishimura,  Proc.  Natl.  Acad.  Sci.  U.S.  I9  (I960) . 

®A.  Terenin,  B.  Putzelko,  I.  Akimov,  Discussions  Faraday  Soc.  83 

(1959). 
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of  cytochrome  do  not  overlap.'^  However,  there  are  no  data  at  present 
which  indicate  which  of  the  three  explanations  given  above  Is  closest 
to  the  truth. 

Much  work  has  gone  into  the  charge  and  energy  transfer  problem  and 
into  Investigations  of  the  semiconductor  properties  of  proteins.®'®'^® 

It  has  been  found  possible  to  calculate  the  width  and  separation  of 
valence  and  conduction  bands  for  certain  configurations  of  proteins.® 
Experiments  on  a  number  of  proteins  i)>  the  polycrystalline  or  dry  gel 
form  have  shown  that  the  electrical  conductivity  is  proportional  to 
exp(-2'E/kT)  where  AS  is  an  activation  energy  of  2.6  to  3-1  ev.^®  How¬ 
ever,  it  is  not  clear  whether  the  band  model  of  semiconductivity  is  valid 
for  systems  such  as  protein  or  amino  acid  molecular  crystals.  Nor  is  it 
clear  that  this  model  is  pertinent  to  charge  transport  within  protein 
molecules.  One  can  argue  against  the  suitability  of  the  band  model  for 
proteins  and  amino  acid  crystals  in  the  following  way.  Molecules  such  as 
amino  acids,  proteins  and  the  like  do  not  form  ionic  or  covalent  crystals 
but  rather  amorphous  solids  or  crystals  in  which  the  molecules  are  held 
to  one  another  by  weak  forces  such  as  the  van  der  Waals  force .  There  is 
but  slight  interpenetration  of  molecular  charge  clouds  at  the  van  der  Waals 
radius.  As  a  result,  any  Intermolecular  conduction  band  capable  of 
supporting  charge  transfer  between  molecules  will  be  very  narrow.  Even 
when  the  molecules  are  linked  by  hydrogen  bonds, where  the  electron 


^Th.  Forster,  Discussions  Faraday  Soc.  7  (1959) • 

®A.  Szent-Gyorgyi,  Nature  lj+8,  157  (19^1)* 

®M.  G.  Evans  and  J.  Gergely,  Biochim.  Biophys .  Acta  188  (1949)- 
^°D.  Bley  and  D.  Spivey,  Trans.  Faraday  Soc.  l432  (I96O) . 

^*■1.  Orgel,  Rev.  Mod.  Phys.  1^0  (l959)  • 
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overlap  Is  enhanced  by  the  presence  of  the  intermediate  proton,  the 
Interaction,  hence  the  band  width,  is  very  small. 

The  velocity  of  an  electron  in  a  conduction  band  is  the  group  velocity 
of  a  wave  packet,  which  is  (l/fi)  (3E^/9k)  ,  where  is  the  energy  as  a 

function  of  the  wave  vector  k  .  Since  the  range  of  k  in  each  band  is 
2jt/a  (a  ■  lattice  parameter)  the  narrower  the  band,  the  smaller  is  the 
average  value  of  9E  /5k  and  the  smaller  is  the  velocity.  In  order  for 
an  electron  to  "see"  the  periodicity  of  the  lattice,  it  must  traverse 
several  unit  cells  between  phonon  collisions.  If  the  electron  velocity 
is  so  low  that  this  does  not  happen,  the  Bloch  functions  are  no  longer 
good  zero-order  approximations  for  the  electron-wave  fiinctions  and 
ordinary  conductivity,  i.e.,  band  theory  is  inappropriate.  In  this  case 
it  seems  better  to  invert  the  usual  orders  of  approximation,  treating  the 
electron-phonon  interaction  in  zero  order  to  get  a  self-trapped  charge 
and  then  Introducing  the  periodic  potential  of  the  crystal  or  molecule 
as  a  perturbation  which  gives  rise  to  a  random  walk  migration  of  the 
trapped  charge. 

In  a  recent  review  article^®  Garrett  points  out  that  the  observed 
temperature  dependence  of  conductivity  in  organic  crystals  may  be  inherent 
in  the  charge  mobility,  instead  of  the  number  of  charge  carriers .  The 
random  walk  migration  of  trapped  charges  does  give  rise  to  a  temperature 
dependence  of  mobility  of  the  form  exp(-AE/Kr)  where  AE  is  a  thennal 
activation  energy.  Garrett  also  mentions  that  no  Hall  effect  has  been 
observed  in  organic  semiconductors  (as  of  1959)  and  that  Seebeck  effect 
investigations  indicate  the  charge  carriers  are  positive. 


^^C.G.B.  Garrett  in  Semiconductors  (edited  by  N.  B.  Hannay,  Relnhold, 
New  York,  1959) • 


-  9  - 


loffe^^  has  suggested  that  the  observed  conductivity  in  organic  solids 
can  be  described  in  termsof  the  mobility  of  ions  in  solid  electrolytes. 

Such  a  description  gives  rise  to  the  observed  temperature  dependence, 

Charge  conduction  within  a  polymer  molecule  such  as  the  a-helix 
polypeptide  (see  Fig.  2),  if  it  occurs,  is  probably  confined  either  to  the 
polymer  back  bone  or  to  the  hydrogen  bond  system.  In  either  case  the 
conducting  charge  must  follow  a  crooked  path,  being  scattered  into  a  new 
direction  after  traveling  only  a  few  Angstroms.  The  momentum  vector  k 
Is  not  a  good  quantum  ntmiber  for  charges  suffering  such  scattering.  Their 
behavior  is  best  described  by  localized  wave  functions  which  migrate  along 
the  polymer.  Thus  the  band  model  also  seems  inappropriate  for  the 
description  of  intramolecular  charge  transport  in  proteins. 

Additional  information  and  references  on  charge  and  energy  migration 
in  biological  systems  may  be  found  in  volume  QJ  of  the  Discussions  of 
the  Faraday  Society. 

It  has  been  known  for  sometime  that  irradiation  with  ultraviolet 
light  causes  electronic  excitations  in  biological  molecules.  The  events 
which  take  place  in  these  molecules  subsequent  to  irradiation  often 
Involve  energy  or  charge  transfer.  When  the  ultraviolet  irradiation  is 
performed  at  low  temperature,  such  as  77°  K,  these  events  are  slowed  or 
modified  and  may  be  observed  experimentally.  Thus  observations  such  as 
the  thermoluminescence  and  electron  paramagnetic  resonance  experiments 
reported  in  this  thesis  shed  some  light  on  charge  and  energy  migration 
processes . 


F.  Ioffe,  Physics  of  Semiconductors  (infosearch,  London,  i960). 
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Thermoluminescence  Is  the  name  given  to  visible  and  ultraviolet 
light  released  from  a  substance  as  a  result  of  heating.  This  phenomenon 
has  been  known  for  about  3OO  years .  Vincenzo  Cascarlolo,  a  cobbler  In 
Bologna,  may  have  observed  It  as  early  as  Cascarlolo,  an  alchemist 

by  avocation,  found  some  particularly  Interesting  rocks  which  sparkled 
in  the  sunlight  and  proceeded  to  heat  them  in  hie  furnace.  In  doing  so, 
he  apparently  observed  the  thermal  release  of  phosphorescence  instead  of 
the  formation  of  precious  metals  as  he  had  hoped.  There  are  reliable 
reports  that  Robert  Boyle  observed  thermoluminescence  from  heated  diamonds 
in  1663.  And  in  I681,  the  Berlin  physician  Johann  Elsholtz  published  an 
account  of  experiments  in  which  fluorspar  powder  was  formed  into  varlovis 
symbols  on  a  metal  plate,  then  heated  and  made  to  glow. 

These  investigators  noticed  that  the  substances  did  not  glow  when 
cooled  and  reheated.  They  realized  that  heating  released  light  from  the 
substance  but  that  heat  was  not  transformed  into  light. 

For  sometime  it  was  believed  that  the  heating  of  thermoluminescent 
materials  released  trapped  sunlight.  However,  it  eventtially  became  well 
known  that  the  color  of  the  luminescence  was  a  characteristic  of  the 
material  and  that  no  luminescence  appeared  unless  the  wavelength  of  the 
exciting  light  was  somewhat  shorter  than  the  characteristic  l\Minescence 
wavelength  for  that  material.  Toward  the  end  of  the  l8th  Century  it  was 
reported  that  various  materials  would  exhibit  thermoluminescence  after 
exposure  to  cathode,  radium  and  other  rays  as  well  as  to  visible  light. 
This  information  gave  rise  to  the  idea,  essentially  correct,  that  the 
radiation  produced  some  kind  of  change  or  excitation  in  the  material, 


^■^R.  W.  Wood,  "Fluorescence  and  Phosphorescence"  in  Encyclopedia 
Britannlca . 
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that  heating  pennitted  the  material  to  relax  or  change  hack  to  the  normal 
etate  and  that  light  was  given  off  during  this  relaxation. 

In  this  century,  the  modern  band  theories  of  solids  have  permitted  a 
q,uantltative  understanding  of  thermoluminescence  from  inorganic  solids, 
particularly  as  a  result  of  the  experimental  work  of  Urbach  and  his 
collaborators^®  in  1930  and  by  Randall  and  Wilkins^®  during  the  early 
19^0 's.  This  work  has  been  summarized  in  some  detail  by  Wood^*  and  by 
Garlick.^’^ 

The  study  of  thermoluminescence  from  organic  or  biological  systems 
appears  to  have  a  much  shorter  history.  Arnold  and  Sherwood  observed 
themoluminescence  from  chloroplasts  which  had  been  irradiated  with  white 
light  at  room  temperatixre . ^ ®  Chloroplasts  are  subcellular  bodies  which 
contain  chlorophyll,  cytochromes  and  a  number  of  enzymes  important  for 
photosynthesis.  Arnold  and  Sherwood  thus  demonstrated  that  chloroplasts 
could  store  energy  received  in  the  form  of  visible  light. 

Augenstine  and  co-workers^®'' studied  thermoluminescence  from  amino 
acids  and  proteins  in  an  attempt  to  understand  the  effects  of  x-ray 
irradiation  on  these  compounds.  They  cooled  powdered  crystalline  samples 
to  liquid  nitrogen  temperature  and  irradiated  them  with  Cobalt-60  gamma 


^®P.  Urbach,  ^  Sitzung  Berichte,  Vienna  Academy  of  Science  (l930) 

^®J.  T.  Randall  and  M.  H.  Wilkins,  Proc.  Roy.  Soc.  A  l84,  366  (19^5). 

Garlick,  "Luiminescence"  in  Encyclopedia  of  Physics  (edited 
by  S.  Flxigge,  Volume  XXVI,  1-128,  1958). 

^®W.  Arnold  and  H.  K.  Sherwood,  Proc.  Natl.  Acad.  Sci.  4-3 1  105  (l957)j 
J.  Phys.  Chera.  7  (1959). 

G.  Augenstine,  J.  G.  Carter,  D.  R.  Nelson  and  H.  P.  Yockey, 
Radiation  Research  Supplement  2,  19  (i960). 

®°C.  J.  Weinberg,  D.  R.  Nelson,  J.  G.  Carter  and  L.  G.  Augenstine, 

J.  Chem.  Phys.  2869  (1962). 
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rays.  After  irradiation,  the  samples  were  heated  at  a  linear  rate  and  the 
thermolminesoence  intensity  as  a  function  of  time  or  temperature  was 
measured.  This  function  is  known  as  the  glow  curve. 

Augenstine  et  have  obtained  the  following  results,  (l)  For 

amino  acids  containing  ring  structures  thermal  activation  energies  for  the 
various  peaks  in  the  glow  curve  are  near  0.2  ev  while  for  proteins,  values 
are  less  than  0.1  ev.  (2)  Thermolumlnescence  intensity  from  amino  acids 
which  contain  ring  structiore  is  as  much  as  three  orders  of  magnitude 
greater  than  that  from  nonring  amino  acids.  (3)  In  amino  acid  crystals, 
the  chemical  structure  is  more  closely  correlated  with  the  characteristics 
of  the  thermolumlnescence  than  is  the  intermolecular  arrangement  as 
reflected  in  such  a  parameter  as  the  crystal  space  group. 

The  thermolumlnescence  work  reported  in  this  thesis  had  two  goals: 

(l)  a  theory  of  thermolumlnescence  applicable  to  the  materials  studied, 
and  (2)  knowledge  of  the  behavior  of  the  ultraviolet  excited  charge  from 
a  study  of  the  glow-curve  kinetics.  A  development  of  kinetics  and  a 
discussion  of  possible  mechanisms  for  themolumlnescence  are  contained 
in  Chapter  II.  Results  of  the  thermolimilnescence  experiments  are  presented 
in  Chapter  III.  Chapter  IV  contains  a  sketch  of  the  theory  of  electron 
paramagnetic  resonance  together  with  the  results  of  the  experiments. 
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CHAPTER  II 


THEORY  AND  KINETICS  OP  THERM0LUMINE3CENCB 

In  this  chapter  three  theoretical  expressions  for  thermolvunineseenee 
intensity  will  be  developed.  The  distinguishing  characteristic  of  each  of 
these  developments  is  the  choice  of  the  physical  model  of  the  luminescent 
material  on  which  the  development  is  based.  The  model  usually  used^’® 
and  the  one  to  be  described  first  is  the  crystal  model,  involving  valence 
and  conduction  bands,  plus  various  electron  and  hole  trapping  sites.  Next 
to  be  developed  is  a  simple  molecular  model  in  which  each  molecule  absorbs, 
stores  and  emits  energy  independently  of  its  environment .  This  model 
proves  useful  in  interpreting  thermoluminescence  of  aromatic  amino  acids. 

The  third  model  Is  a  polymer  or  one-dimensional  crystal  model  in  which 
trapped  charges  may  hop  from  one  monomer  or  unit  cell  to  an  adjacent  one. 

This  third  point  of  view  is  proposed  to  describe  thermoluminescence  from 
polypeptides  and  proteins.  In  each  case,  the  nature  of  the  energy  trapping 
or  metastable  excited  state  is  discussed  but  not  specified  precisely  and 
the  form  of  the  temperature-dependent  probability  for  thermal  depopulation 
of  these  states  is  assumed.  It  is  also  assvaned  that  radiative  decay  to 
the  groiind  state  follows  thermal  emptying  of  a  metastable  excited  state 
with  10  sec .  The  kinetics  of  the  luminescence  is  then  derived  on  the 
basis  of  the  model  and  these  Initial  assumptions. 

A  theory  of  thermoluminescence  from  crystalline  materials  is  given 
by  Ai;igenstine,  ^  They  present  data  on  the  x-ray- induced  glow  curves 

^See  for  example,  D.  Curie,  Luminescence  Crlstalline  (Dunod,  Paris,  I960) . 

®G.F.J.  Garlick,  "Lxminescence"  in  Encyclopedia  of  Physics  (edited 
by  S.  Flugge,  Volume  XXVI,  1-128,  1958). 

^L.  G.  Augenstine,  J.  G.  Carter,  D.  R.  Nelson,  and  H.  P.  Yockey, 

Radiation  Research  Supplement  2,  I9  (1960) . 
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of  proteins  and  amino  acids  Interpreted  in  terms  of  their  theory.  We  wish 
to  give  a  more  detailed  theory,  based  to  some  extent  on  the  work  of 
Halperln  and  Braner.*  Their  notation  will  be  used  here. 

To  begin  with,  we  assume  that  our  sample  Is  a  crystalline  material 
in  which  the  electrons  are  subjected  to  a  periodic  potential  so  that  the 
Bloch  factions  are  appropriate  as  electron-wave  functions  and  the  energy 
levels  of  the  electrons  separate  into  bands.  Also,  we  suppose  that  there 
is  a  gap  between  the  top  of  the  highest  band  of  filled  states,  the  valence 
band,  and  the  lowest  band  of  empty  states,  the  conduction  band.  We  assume 
that  the  width  of  this  gap  is  much  greater  than  kT  at  room  temperature. 
Because  of  the  existence  of  lattice  imperfections,  randomly  distributed 
dipoles,  etc.,  there  will  be  localized  states  available  to  the  electrons 
and  holes  within  the  energy  gap.  The  energy  E  is  measured  from  any 
convenient  level,  such  as  the  top  of  the  valence  band.  The  magnitude  of 
the  energy  gap  will  then  be  E  ,  as  shown  in  Pig.  4. 

After  ultraviolet  or  x-ray  irradiation,  a  number  of  electrons  will 
have  been  excited  into  the  conduction  band  and  subsequently  trapped  in 
the  localized  electron  states,  called  traps  from  now  on.  The  holes  left 
behind  in  the  valence  band  will  have  migrated  to  the  localized  hole  states, 
called  limiinescence  centers.  There  are,  say,  n^(E)dE  trapped  electrons 
per  unit  volume  at  energy  E  to  S  +  dE  ,  m(E’)dE’  holes  per  unit 
volume  at  energy  E'  to  dE'  .  As  the  sample  is  heated,  the  trapped 
electrons  and/or  holes  acquire  sufficient  thermal  energy  to  escape  into 
the  conduction  or  valence  band.  Once  excited  into  bands,  the  charges 
migrate,  and  the  electrons  and  holes  recombine.  In  describing  this  process, 

^A.  Halperln  and  A.  A.  Braner,  Phys.  Rev.  117 >  408  (i960). 
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FIG.  4 --Energy  level  diagram  indicating  valence  and  conduction 
bands  of  states  plus  the  relatively  few  states  or 
traps  In  the  energy  gap. 
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we  assume  that  the  electrons  are  the  thermally  activated  particles,  the 
holes  Training  trapped.  The  kinetics  are  the  same  if  the  roles  are 
reversed. 

Before  deriving  these  kinetics,  we  consider  the  equilibria  distri¬ 
bution  of  electrons  in  the  crystal.  Let  N(b)  be  the  density  of  electron 
states  throughout  the  valence  band,  conduction  band,  and  energy  gap.  Let 
f(E)  stand  for  the  Fertni-Dlrac  distribution 

"  exp[(E  -  E^/kT]  +  1  ' 


in  which  E^  Is  the  Fermi  energy. 

Let  n(E)  =  N(E)f(E)  stand  for  the  equilibrium  distribution  of 
electrons  in  the  crystal.  Next,  assume  the  energy  E^  is  located  within 
the  energy  gap  sufficiently  far  from  any  trapping  levels  so  that,  to  a 
good  approximation, 

f(E)=expf-^ - (2-2) 

\  KT  / 


n(E)  =  N(E)exp|^- 


KT  / 


(2-3) 


Consider  two  thin  constant  energy  shells  of  electron  states,  one 

about  E.  in  the  energy  gap  and  the  other  at  B  near  the  bottom  of 
X  c 

the  conduction  band  as  shown  in  Fig.  5*  The  equilibrium  electron  densities 
in  these  two  shells  are 


n(E^)  =  N(E^) 
n(E^)  =  N(E^)  exp 


(2-4) 
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PIG,  5~*Energy  level  diagram  indicating  the  two  thin  shells 

of  states,  one  in  the  energy  gaps,  one  in  the  conduction 
band,  as  described  in  the  text. 


N(B^)  -  ■  number  per  unit  volume  of  unoccupied  trap  states  In  the 

interval  to  B^  +  dE^  , 

N(B^)  -  n(B^)  ■  number  per  unit  volume  of  unoccupied  conduction  states 

in  the  interval  E  to  E  +  dE  , 
c  c  c 

The  principles  of  detailed  balance  tells  us  that  at  equilibrium  the 

number  of  electrons  excited  from  the  shell  at  E.  to  the  shell  at  E 

x  c 

per  unit  time  is  equal  to  the  number  trapped  from  B  to  E,  per  unit 

C  w 

time.  This  fact  can  be  lased  to  find  the  temperature  dependence  of  the 

rate  of  excitation  of  electrons  from  E.  to  E  .  Let 

X  c 

73^(t)[N(eJ  -  n(E^)] 

be  the  probability  per  unit  time  that  an  electron  in  a  trap  at  E.  will 

X 

be  excited  to  the  conduction  band  at  E  .  Let 

c 

-  n(E^)] 

be  the  probability  per  unit  time  that  an  electron  in  the  conduction  band 
will  fall  into  one  of  the  [N(E^)  -  n(E^)]  empty  traps  available  in  a 
imlt  volume.  We  assume  is  not  dependent  on  temperature.  The 

parameters  ^^^(t)  ,  7g  contain  such  quantities  as  the  matrix  element 
for  the  transition  Involved  and  the  energy  distribution  of  phonons  which 
are  involved  in  the  transition.  Attempts  at  a  calculation  of  these 
quantities  have  been  made  by  several  authors.^'® 

Prom  the  statement  that  the  number  of  electrons  excited  per  unit 
time  equals  the  number  trapped  we  get  the  equation 

n(E^)7i(T)[N(E^)  -  n(E^)]  =  n(B^)7g[Il(E^)  -  n(B^)]  .  (2-5) 


®R.  Kubo,  Phys.  Rev.  86,  929  (1952). 

®B.  Goodman,  A.  W.  Lawson,  L.  I.  Schiff,  Phys.  Rev.  Jl,  174  (194t). 
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Using  Eq.  (2-4)  and  cancelling  like  factors  from  both  sides  we  can 
rewrite  Bq.  (2-5)  as 


yi(T) 


or,  consistent  with  approximation  (2-2),  as 


^2.(1)  =  72 


On  the  basis  of  this  result  we  will  take 


7g  =  constant  =  y 


7j^(t)  =  7  expl- 


E  -  E^' 

_c _ t 

KT  , 


(2-6) 


and  assume  these  results  are  valid  for  any  distribution  of  electrons  in 
the  crystal.  In  the  appendix  we  outline  a  calculation  of  7^(t)  due  to 
Kubo.®  Using  the  Einstein  model  of  a  crystal  and  certain  additional 
assumptions,  he  obtains  the  same  functional  form  of  7^(T)  as  (2-26) 
and  an  expression  for  y  . 

The  simplest  glow  c\rrve  has  one  peak  and  all  of  the  emitted  photons 
have  the  same  energy.  Such  a  curve  can  be  described  on  the  basis  of  one 
trapping  level  at  energy  E^  ,  traps  per  unit  volume,  and  one  level 
of  trapped  holes.  At  the  end  of  the  irradiation  period,  let  us  say  there 
are  n  trapped  electrons  per  unit  volume.  Since  the  probability  for 
thermal  excitation  of  a  trapped  electron  to  energy  E'  within  the 
conduction  band  is  proportional  to 


exp 


; 


nearly  all  thermally  excited  electrons  will  go  to  states  near  the  bottom 
of  the  conduction  band  where  S'  =  E 

c 
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It  is  helpful  to  modify  notation  at  this  point  to  discuss  changes 


/ 


n 


In^l 


in  electron  and  hole  concentration  with  time.  Let  n^  «=>  n(B^)  ,  ' 

n^  =  follows,  t  refers  to  traps  when  it 

appears  as  a  subscript  and  to  time  when  it  appears  as  the  argument  of  a 


function.  Thus  n,  (t)  is  the  density  of  electrons  . in  traps  at  energy 
B,  at  time  t  .  In  this  notation,  the  change  in  the  concentration  of 


trapped  electrons  during  time  dt  is 

/  E. 


dn,  =  -  n  7(N  -  n^)  exp  - 


KT 


dt  +  n^7(N^  -  n^)dt  .  (2-7) 


The  change  in  the  concentration  of  conduction  electrons  during  this 
interval  is 

dn  =  -  dn.  -  n  Arndt  , 
c  t  c  ' 

where  A  is  a  constant  depending  on  the  recombination  efficiency  and 
m(t)  is  the  number  of  holes  per  unit  volume.  The  change  in  the  concen¬ 
tration  of  holes  during  this  interval  is 

dm  =  -  n  Arndt  . 
c 

The  set  of  differential  equations  describing  the  kinetics  of  return  to 
equilibrium  distribution  is 


dn  /  E  -  E  A 

~  =  -  n^r(N^  -  n^)  expU 


dn^  dn^  dm 

dt  dt  dt 


dm  , 

dt  =  -  “^t  • 


(2-8) 

(2-9) 

(2-10) 


Braner  and  Halperin*  give  a  careful  discussion  of  the  solution  to 
(2-8),  (2-9)  and  (2-10)  using  the  approximation  dn^/dt  =  0  . 
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P.  K.  Weyl"^  derives  kinetics  for  thermoluminescence  using  a  different 
approximation.  He  replaces  n^(t)  with  what  he  calls  the  thenaodynamio 
eq,uilibrium  value 

/  E  - 

n.  (t)  expl-  - 

V  kT 


This  replacement  is  Incorrect  since  it  does  not  take  into  account  the 
difference  in  the  density  of  states  between  the  energy  gap  and  the 
conduction  band. 

The  approximation  we  have  developed  begins  by  assuming  first  that 
the  states  of  the  conduction  band  are  always  largely  \moccupied,  second 
that  were  it  not  for  recombination  of  conduction  band  electrons  with 
holes,  n  and  n  would  be  in  approximate  thennal  equilibrium.  The 

C  U 

relative  rates  of  recombination  and  retrapping  depend  upon  the  ratio 
A/j  .  The  greater  this  ratio,  the  more  will  n^  be  reduced  from  its 
equilibrium  value  with  respect  to  n^  . 

On  the  basis  of  this  reasoning  we  put 


in  the  second  term  on  the  right-hand  side  of  Eq. 
approximations  change  (2-8)  to 

1  -  exp  ( 


dn  /  E  > 

-  =  -  n  7N  expl - 

dt  ^  \  kT, 


(2-8). 


(2-11) 


These  two 


p  yN  /AN  /  E  \ 

+  n^  — ^exp(--j  expf - j  (2-12) 

"t  ^ 


where  E  =  E  -  E  from  now  on-  This  equation  is  of  the  form  known  as 
c  X 


■^P.  K.  Weyl,  J.  Chem.  Phys.  26,  (1957). 
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Bernoulli's  equation.  The  solution  Is 


n.(0)  exp[-/Q  B(t')dt'] 

n  (t) - ^ - 

1  -  n^(0)  /q  C(t' )  ejcpC-  /J  B(t")dt"]dt' 

B(t)  »  7N^[i  -  exp(-  |Jexp[-  3^^ 

yN  /  A\  r  E  ■ 

C(t)  =  — expi-  —  j  exp  -  — — 

N.  V  7/  L  l5;T(t)_ 


(2-13) 


(2-14) 


(2-15) 


Equation  (2-9)  can  be  solved  formally,  giving 

n^(t)  =  m(t)  -  n^(t)  +  m(o)(p  -  l) 


(2-16) 


where 


n^(0)  +  n^(0) 


(2-17) 


and  n  (0)  ,  n  (O)  ,  m(0)  are  the  initial  electron  and  hole  concentrations, 

w  C 

Substituting  (2-I6)  into  (2-10)  we  have 


=  m(t)  A[n.  (t)  -  (p  -  l)m(o)]  -  Am^(t)  . 


(2-18) 


Equation  (2-l8)  is  also  a  Bernoulli  equation.  Its  solution  is 

m(0) 

m(t)  =  - i— -rvv\ - 


1  +  Am(0)  /q  e^^^  '  dt' 


(2-19) 


where 


u(t)  ^  ^  J'  [*^(t')  -  ni(o)(p  -  l)]dt' 


(2-20) 


and  n^(t)  is  given  by  (2-I3). 

The  Intensity  of  the  thermolurainescence  will  be  proportional  to  the 
rate  of  recombination,  i.e.,  to  dm/dt  .  Letting  g  be  the  proportion¬ 
ality  constant  which  depends  upon  the  geometry  of  the  experiment,  the 
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detection  efficiency  of  the  photomultiplier  tube,  etc.,  we  have 


I(t) 


gAm(O)  e^^^^ 

1  +  Aii(O)  /q  e’^^^'^dt* 


X 


n^(t)  +  (l  -  p)in(o) 


,(0) 

1  +  Jaaio)  /q 


.(2-21) 


Expression  (2-2l)  gives  the  luminescence  from  a  crystalline  sample  which 
has  at  time  zero  h^(0)  trapped  electrons  and  m(o)  recombination  centers 
per  unit  volume  as  a  result  of  some  previous  irradiation. 

We  turn  now  to  a  model  for  thermoluminescence  based  on  a  homogeneous 
sample  of  weakly  interacting  molecules.  In  this  model,  each  molecule  is 
presumed  to  act  independently  in  the  processes  of  absorbing,  storing  and 
re-emitting  energy  in  the  form  of  a  photon. 

There  are  several  ways  in  which  a  molecule  could  trap  and  store 
energy.  Perhaps  the  most  familar  is  the  radiationless  conversion  from  an 
excited  singlet  to  excited  triplet  state  which  occurs  in  many  aromatic 
molecules.  The  triplet  state  is  metastable  since  the  spin  of  the  excited 
electron  must  be  flipped  during  the  triplet-singlet  transition  and  the 
spin-orbit  interaction  which  promotes  this  flipping  is  weak  in  these 
molecules.  There  are  two  difficulties  with  the  hypothesis  that  a  molecular 
triplet  state  is  the  metastable  or  trapping  state  required  for  thermo- 
luminescence.  One  is  the  shortness  of  lifetime.  Calculations  of  the 
lifetime  of  the  triplet  state  in  aromatic  molecules  have  been  made  by 
McClure.®  Basing  the  calculation  on  the  estimated  spin-orbit  coupling 
provided  by  the  atom  nearest  the  orbitals  Involved  in  the  transition  he 


®D.  S.  McClure,  J.  Chem.  Phys.  1^,  905  (19^9)/  20,  682  (l952). 


obtains  lifetimes  of  0.003,  0,0004,  0.0001  seconds  when  the  coupling  is 
due  to  a  carbon,  nitrogen  or  oxygen  atom  respectively.  When  molecular 
symmetry  is  taken  in  account  there  is  seme  cancellation  of  this  coupling 
and  an  increase  in  lifetime  to  as  much  as  several  seconds.  The  phenomenon 
of  thermoluminescence,  however,  requires  a  metastable  state  with  a  life¬ 
time  several  orders  of  magnitude  longer,  particularly  at  low  temperature. 
This  is  the  second  difficulty.  There  is  no  theoretical  Justification  for 
an  assumption  that  the  lifetime  of  the  triplet  state  is  strongly  tempera¬ 
ture  dependent.  Thus  the  triplet  state  does  not  have  the  necessary 
properties  with  regard  to  lifetime  to  be  considered  the  metastable  excited 
state  involved  in  themolimilnescence. 

One  possibility  for  a  molecular  excited  state  with  the  required 
properties  is  the  trapping  of  an  excited  electron  by  an  ionic  region  of 
the  molecule.  For  example,  amino  acids  in  water  solution  at  physiological 
pH  have  both  the  amino  and  carboxyl  groups  ionized.  Some  small  percentage 
of  the  molecules  may  also  have  this  configuration  in  the  dry  state  and 
thus  have  the  ability  to  form  the  necessary  metastable  excited  state. 

Let  us  call  this  conjecture  the  ionic  trapping  hypothesis. 

Another  possibility  is  suggested  by  the  properties  of  certain  ethylenlc 
compounds.®  These  molecules  can  form  cls-trans  isomeric  states  as  the 
result  of  rotation  about  a  carbon-carbon  double  bond.  The  potential  energy 
of  the  two  electrons  involved  in  the  «  orbital  of  this  double  bond  is 
a  function  of  the  angle  of  bond  rotation  as  shown  in  Fig.  6.  The  curve 
labeled  S  gives  the  energy  of  the  singlet  state;  the  curve  labeled  T 
gives  the  energy  of  the  triplet  state  as  a  function  of  bond  rotation  angle. 


®C.  Reid,  Excited  States  in  Chemistry  and  Biology  (Butterworth, 
London,  1957)- 
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Angle  about  double  bond 


FIG.  6--Energy  of  singlet  (s)  and  triplet  (T)  state 
orbitals  as  a  function  of  the  angle  of  twist 
about  the  axis  of  the  ethylenlc  double  bond 
(after  Reid®) . 


The  molecules  Investigated  are  in  the  singlet  state  before  and  after 
isomerization.  Some  but  not  all  of  them  appeeu*  to  undergo  a  singlet- 
triplet-singlet  transition  during  the  isomerization. 

The  reaction  rate  for  the  isomerization  has  the  form 

K  =  A  exp(-  1^)  .  (2-22) 

On  the  basis  of  experiments  which  determine  A  and  E  ,  the  molecules 
investigated  may  be  divided  into  two  groups,  (l)  A  group  for  which  the 
activation  energy  is  low  and  the  frequency  factor  A  is  abnormally  low, 
being  about  10  sec  .  (2)  A  group  for  which  the  activation  energy  is 
about  twice  as  large  as  in  the  first  group  and  the  frequency  factor  is 
much  closer  to  values  found  in  most  reactions,  l.e.,  about  10^^  .  The 
interpretation  is  that  the  first  group  of  molecules  undergoes  the 
isomerization  via  the  singlet -triplet -singlet  pathway.  The  required 
activation  energy  for  this  pathway  is  relatively  low,  and  the  forbidden¬ 
ness  of  the  singlet-triplet  transition  accounts  for  the  low  frequency 
factor.  The  second  group  isomerizes  by  going  over  the  singlet  state 
energy  barrier. 

An  interesting  feature  of  the  first  group  is  that  at  sufficiently 
low  temperatures,  molecules  may  become  trapped  in  the  intermediate  90° 
state  at  position  C  .  The  probability  for  sucii  molecules  to  become 
untrapped  is  proportional  to 


and  thus  has  the  temperature  dependence  needed  for  therraoluminescence. 
If  such  an  isomerization  effect  is  to  explain  themnolumlnescence  from 
amino  acids,  two  additional  assumptions  must  be  made.  One  is  that  an 
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analogous  isomerization  effect  may  occur  in  some  amino  aeldS;  say  at  the 
carbon-carbon  bond  linking  the  side  group  to  the  Of-carbon.  The  second  is 
that  a  photon  will  be  given  off  at  some  point  as  the  molecule  goes  from 
the  metastable  state  (at  position  C)  to  the  ground  state  at  0°  to  l80°. 
There  is  no  evidence  at  present  to  substantiate  either  assumption.  We 
call  this  possible  explanation  for  amino  acid  thermolumineseence  the 
frozen  isomerization  hypothesis. 

Thirdly,  it  is  possible  that  the  necessary  molecular  trapping  state 
is  provided  by  a  positive  metal  ion  which  is  complexed  with  the  molecule. 

A  number  of  metal  ions,  e.g.,  iron,  copper,  have  a  strong  tendency  to 
complex  with  amino  acids.  Irradiation  of  such  a  complex  at  liquid 
nitrogen  temperature  may  cause  the  formation  of  a  stable  charge-transfer 
complex.  When  the  temperatvire  is  raised  the  transfered  electron  reverts 
back  to  an  excited  state  of  the  molecule.  From  this  state  it  undergoes 
a  transition  to  the  ground  state  with  the  emission  of  a  photon.  Let  us 
call  this  the  charge -transfer-complex  hypothesis. 

To  derive  an  expression  for  the  thermolumineseence  intensity  on  the 
basis  of  a  molecular  model  such  as  discussed  above  we  begin  by  postulating 
that  the  transition  rate  for  a  molecule  going  from  a  metastable  excited 
state  to  a  short-lived  (singlet  or  triplet)  excited  state  and  then  to 
the  ground  state  with  the  emission  of  a  photon  is 

w  =  7  exp(J-  (2-23) 


where  E  is  an  activation  energy  and  where  7  contains  the  matrix 
element  for  the  electronic  transition.  If  there  are  n  molecules  per 
unit  volume  in  the  excited  state,  the  number  of  radiative  transitions 
expected  during  time  dt  is 


-  dn  =  n7  exp 


28 


dt  . 


(2-24) 


The  solution  to  this  differential  equation  Is 


n(t)  »  n(0)  exp  ^ -  7 


7  f  e*p[-  5^]  at 'I 


(2-25) 


The  intensity  of  the  luminescence  will  be  proportional  to  dn/dt  .  Thus 

I(t)  s  gn(0)7  exp^-  ^  j  exp -jj  y  J  expj^-  •  (2-26) 

The  parameters  E  ,  y  which  appear  in  this  expression  can  be  determined 
from  the  experimental  curves  In  the  following  way.  Let  time  t  =  0 
correspond  to  the  beginning  of  the  themoluminescence .  Then  for  times 
corresponding  to  the  initial  rise  of  the  glow  curve,  the  integral  in 
(2-26)  will  be  small  and  the  luminescence  will  be  given  approximately  by 


r(t)  =  gn(0)7  exp^- 


(2-2T) 


Using  this  region  of  the  glow  curve  for  a  plot  of  In  l(t)  vs  [l/T(t)] 
one  obtains  a  straight  line  with  slope  =  -  E/k  .  At  the  peak  of  the 


glow  curve, 


T  =  T 


^  dt  l,dt/ 


and  4^  =  0  . 
dt 


Differentiating  (2-26)  and  introducing  these  peak  conditions,  one  finds 


*  fdT\  f  E 

^  “  .  ^2  U/P  \ 


(2-28) 


These  results  will  be  used  in  the  analysis  of  the  experiments  discussed 
in  the  next  chapter.  Equation  (2-28)  reveals  another  feature  of  the 
glow  curves,  namely,  that  the  temperature  at  which  the  peak  occurs  depends 
on  the  heating  rate.  This  is  clearly  seen  in  Fig.  ?•  A  number  of  workers^ 


^D.  Curie,  op.  clt. 


29 


relationship  between  heating  rate  and  peak  temperature. 


have  attempted  to  use  this  fact  to  determine  the  constants  E  ,  y  . 

Unfortunately  this  method  laclts  accuracy  because  T  changes  only 

P 

slightly  with  dT/dt  j  it  is  difficult  to  reduce  the  error  below  20-30?t. 

Now  we  come  to  the  mechanism  proposed  to  account  for  thermolumlnes- 
cence  from  polymers,  more  specifically  from  polypeptides.  In  this  case, 
one  imagines  that  irradiation  of  the  molecules  produces  holes  and  trapped 
charges  within  the  molecules.  Irradiation  causes  a  charge,  say  an  elec¬ 
tron,  to  be  ejected  from  some  site  on  the  peptide j  it  is  trapped  at  some 
other  site.  If  the  vacancy  is  represented  by  an  open  circle  and  the 
trapped  electron  by  a  closed  circle,  the  polypeptide  can  now  be  represented 
as  in  Pig.  8. 

These  traps  may  be  considered  as  regions  of  local  polarization 
within  the  polypeptide.  The  site  of  the  a-carbon  of  the  amino  acid,  is 
likely  to  be  such  a  region.  To  support  this  contention,  there  is  the 
fact  that  amino  acids,  in  water  solution  at  physiological  pH  prefer  the 
ionized  structure  about  the  a-carbon  as  discussed  previously.  Also, 
there  is  work  of  Box,  Freund  and  Lilga^°  showing  that  x-ray  irradiation 
of  some  simple  peptides  produces  dissociation  of  hydrogen  from  the  a-carbon 
atom  of  the  amino  acid,  and  the  work  of  Katayama  and  Gordy^^  showing  that 
the  a-carbon  site  is  a  trap  for  radiation  damage.  Let  us  assvime  that  the 
polar  environment  of  the  a-carbons  in  the  polypeptide  serves  to  trap 
migrating  electrons  or  holes. 

The  transition  of  a  charge  from  one  trap  to  the  next  is  conceived 
of  as  a  barrier  penetration  req,ulring  a  certain  themal  activation  energy. 


^°H.  C.  Box,  H.  G.  Freund  and  K.  Lilga,  Symposium  on  Free  Radicals 
in  Biological  Systems  (ed.  M,  S.  Blois,  Academic  Press,  New  York,  1961) • 

Katayama  and  W.  Gordy,  J.  Chem.  Phys.  117  (1961). 
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FIG.  8--Block  representation  of  a  polypeptide.  Each  square 
represents  a  possible  trapping  site  - 
one  per  amino  acid.  The  black  dot  represents  a 
trapped  electron,  the  open  circle  is  the  corresponding 
hole. 


FIG.  9 — Block  representation  of  polypeptide  showing  the 
ordering  of  trapping  sites  used  in  the  Markov 
analysis  of  charge  migration. 
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It  is  «3gpeeted  that  barriers  within  a  polypeptide  are  more  easily  penetrated 
than  barriers  between  two  different  polypeptides.  For  this  reason,  it  is 
assumed  that  a  charge  will  stay  on  one  molecule,  migrating  along  the 
molecule  until  it  finds  a  hole  with  which  it  may  recombine  to  produce  a 
photon.  It  is  further  assumed,  in  the  Interest  of  simplification,  that 
there  is  only  one  charge  per  molecule  which,  in  fact,  has  been  released 
from  a  trap  site.  Our  model,  then,  consists  of  a  single  electron  and  a 
single  hole  constrained  to  a  linear  array  of  traps  (finite  in  number). 

The  electron  hops  from  one  trap  to  another  but  the  hole  remains  fixed. 

The  process  continues  until  the  electron  and  hole  recombine.  This 
migration  process  has  the  Markov  property:^®  given  the  state  of  the 
molecule  at  any  time,  future  changes  are  independent  of  the  past.  Let 
P^^(T,t)  be  the  conditional  probability  of  finding  the  charge  at  the 
nth  trapping  site  at  time  t  given  that  at  previous  time  the  charge 
was  at  the  ith  trap.  Then 


T  <  .  <  t  .  (8-29) 

m 

Next  we  postulate  that  for  each  trap  in  the  molecule  there  exists  a 
continuous  function  which  is  the  probability  that  an  electron 

could  escape  from  the  rath  trap  during  the  time  interval  t  ,  t  +  6t  . 

Thus 


,,  1-P  (t,t+6t) 

'.So - ^ - • 


(2-30) 


Furthermore,  for  any  two  traps  m  and  n  there  exists  a  function  Pjjjjj(b) 


Gr.  Feller,  Introduction  to  Probability  Theory  and  its  Applications 
(Vol.  1,  2nd  ed.,  Wiley,  New  York,  195T)- 
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which  is  the  conditional  probability  that  if  an  electron  escapes  from 
trap  m  during  t  ,  t  +  8t  it  is  captured  by  trap  n  ,  Thus 

eSo - - <2-3l) 

with 

Zp  =  1  ,  p  “  0  , 
am  '  nmi  ’ 

n 

m  ,  n  finite.  A  stochastic  process  which  conforms  to  these  postulates 
can  be  described  by  the  Kolomogorov  equations 

It  .  yPi.(T,t)C^(t)p.„(t)  .  (S-32) 

In  this  system  of  equations,  1  and  t  are  fixed  and  the  equations  are 
actually  ordinary  differential  equations. 

A  simple  example  will  Indicate  how  the  glow  curve  may  be  obtained 
from  the  above.  Assume  the  molecule  has  one  trap  and  one  recombination 
site,  labelled  1  and  2  respectively.  After  the  material  has  been 
irradiated,  some  molecules  will  have  their  traps  occupied  by  electrons 
(and  recombination  sites  occupied  by  holes) .  The  probability  per  unit 
time  that  an  electron  will  escape  from  the  trap  (when  the  material  is 
warmed)  is  c^(t)  while  the  probability  per  unit  time  that  an  electron 
will  escape  from  a  recombination  site  is  zero,  or  C2(t)  =  0  .  Further, 
Pl2(t)  =  1  ,  that  is  the  electron,  if  it  escapes  from  the  trap,  is 
certain  to  end  up  at  the  recombination  site.  need  not  be  specified. 

Under  these  circumstances,  the  Kolomogorov  equations  take  the  simple 

form 

^  =  C^(t)P^^(T,t) 

It  • 

3^ 


(2-33) 


We  note  that  for  the  8impl.e  case  under  consideration  P^^(T,t)  represents 
•  the  probability  that  a  trapped  electron  has  not  recombined,  Hence  the 

number  of  molecules  per  unit  volume  having  traps  occupied  by  electrons 
at  time  t  is 

n(t)  =  n(T)P^^(T,t)  (2-35) 

where  n(T)  is  the  initial  concentration  of  molecules  with  trapped 
charges.  Hence  Eq.  (2-3^)  may  be  written 

n(t)  =  -  C^(t)n(t)  .  (2-36) 

If  C^(t)  is  known,  Eq.  (2-36)  gives  the  glow  curve  since  the  Intensity 
of  the  emitted  light  is  proportional  to  dn/dt  .  Suppose  that  the 
removal  of  an  electron  from  a  trap  involves  a  thermal  activation  energy 
E  .  Then 

C^(t)  =  y  exp[-  (2-37) 

which  depends  on  time  only  to  the  extent  that  hhe  temperature  T  may 
depend  on  time.  We  have  then,  for  the  glow  cui*ve, 

r  t 

l(t)  =  gn(T)7  exp^-  exp  -  7  j^exp^-  dt’j  (2-38) 

-  T 

which  is  the  same  as  (2-26).  This  is  to  be  expected  since  the  same 
physical  picture  is  used  in  the  derivation  of  both  equations.  For  the 
general  case  we  consider  a  linear  polymer,  £  units  long,  containing 
one  trapped  electron  and  one  recombination  site.  The  electron  may  be 
trapped  in  any  one  of  i-1  sites,  from  which  it  can  escape,  or  in  the 
one  recombination  site  from  which  it  cannot  escape.  Thus,  as  shown  in 
Fig.  9^  the  electron  is  in  trap  i  and  the  recombination  site  is  at  r  . 

•  As  the  material  is  warmed,  it  is  assumed  that  the  electron  may  hop  from 
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1  to  either  i  +  1  or  1  -  1  .  The  process  continues,  by  a  succession 
of  hops  to  neighboring  sites,  until  the  electron  finds  the  site  r  , 
recombination  takes  place  and  a  photon  is  emitted. 

In  order  to  obtain  the  glow  curve  we  need  to  know  dn(t)/dt  where 
n(t)  is  the  total  number  of  molecules  per  unit  volume  containing  a 
trapped  electron  at  time  t  .  We  proceed  as  follows:  Let  be 

the  concentration  of  molecules  at  initial  time  with  electron  at  i  , 
recombination  site  at  r  .  The  number  of  these  molecules  which  at  time 
t  still  have  a  trapped  electron  but  not  necessarily  at  i  ,  is 

-  ni,(T)  Y,  ■  (2-39) 

m<r 

The  number  of  molecules  per  unit  volume,  n(t)  having  an  electron  in  some 
trap  and  a  recombination  site  at  any  of  the  1  possible  sites  is  obtained 
by  summing  over  i  and  r  .  That  is 

L 

n(t)  ^  “^ir^^^  ■  (2-40) 

r*s2  l<r 

The  counting  scheme  used  is  responsible  for  the  summation  over  recombina¬ 
tion  sites  beginning  at  2.  We  assume  that  all  initial  arrangements  of 
trapped  electron  and  recombination  sites  are  equally  probable,  i.e., 

ni^(T)  =  Hq  (2-1^1) 

where  n^  is  a  constant,  the  same  for  all  r  and  all  i<r  . 

Recall  that  the  equations  for  the  conditional  probability  ^^(■’■>t) 

are 

It  =  -  0_(t)P,„(T,t)  +  yPi^(T,t)0^(t)p^„  (2-1.2) 

with  m  *  1,2, ,..,r-l  .  The  assumption  that  when  a  trapped  electron 
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escapes  it  will  be  retrapped  by  one  of  Its  two  nearest  neighbors  gives 

Pjgjj  =  ^  when  |k  -  m  I  o  1  (2-l<-3) 

except 

=  1  (2-44) 

with 

Pto  "  ° 

otherwise . 

Next  we  assume  that  the  escape  probability  Is  the  same  for 

all  traps,  while  for  the  recombination  site  r  ,  C^(t)  =  0  .  Let 

Cjt)  =  C(t)  ,  m  =  l,2,...,r-l  .  (2-46) 

Now  using  (2-43),  (2-44),  (2-45),  (2-46),  we  can  write  the  set  of 
equations  (2-42)  in  vector  form 

-1  i 
1  -1  i 

5  -1 

O 

The  substitution 

Pl^(T,t)  =  U^(t)  exp  -  J  C(t')<it'  (2-48) 

~  T 

simplifies  the  matrix  involved  in  this  equation  and  gives  the  following 
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i 

[ 


n 

:  1 
i,  1 

vector  equation  for  the  U^(t)  .  j  ^ 


0  i 

0l2(t) 

1  0  i  Q 

• 

I-  0  . 

• 

-  c(t) 

«  •  * 

• 

• 

• 

•  •  • 

0  .  .  i 

• 

i  0 

If  we  denote  the  matrix  In  (2-i^9)  hy  ,  this  equation  can  be  written 
more  compactly  as 


U^(t)  =  C(t)B^Ujt) 


The  Initial  conditions  for  the  U.  (t)  are 

im 


u,  (t)  =  6, 

ira'  '  Im 


(2-50) 

(2-51) 


To  solve  (2-50)  we  find  the  matrices  A  ,  A"^  such  that 

r  r 

A"^B  A  =  D  (2-52) 

I*  I*  r  I*  \  / 

where  D  is  the  (r  -  l)  x  (r  -  l)  diagonal  matrix  whose  diagonal 
r 

elements  ^  are  the  eigenvalues  of  B^  .  We  assume  that 

can  always  be  diagonalized.  This  has  been  verified  for  r  S  7  by 
noting  that  ,  the  characteristic  polynomial  of  ,  satisfies  the 
recursion  relation  2  ^  ^  \  ~  ^  > 

and  using  this  relation  to  find  q^  and  its  roots ,  We  now  write  (2-50)  as 

It  A"^U,(t)  =  C(t)A"^B  A  A"V(t) 
otri''  ''rrrrl 
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or 

^  V^(t)  =  C(t)D^V^(t) 

where 

v^(t)  =  A;;\(t)  . 

Equation  (2-53)  can  be  solved  immediately.  We  have 


V^(t)  =  E^(t)V^(T) 


where  E^(t)  is  the  diagonal  matrix 


By(t)  = 


V 


with 


J  =  /  C(t'  )dt' 


The  initial  values  of  the  given  by 


v^(t)  -  a;\(,)  . 

The  solution  to  Eq.  (2-50)  may  now  be  written 


U,(t)  =  A  B  (t)A“V  (t) 
i'  '  r  r'  '  r  1 

and 

PjT,t)  =  e"'^U^(t)  . 

n(t)  may  now  be  computed  from  Eqs.  (2-39)>  (2-^^■0)  and  (2-59)  • 


(2-53) 

(2-5i^) 

(2-55) 


(2-56) 

(2-57) 

(2-58) 

(2-59) 

Specifically, 
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we  have 


n(t)-n(T)^  (2-6o) 

rag  l_i<r  m<r 

where  the  double  sum  Inside  the  parenthesis  is  equal  to  the  sum  of  all 

the  elements  in  A  B  (t)A"^  . 

r  r  r 

The  special  case  L  a  4  provides  an  example  of  this  technique. 
Consider  first  the  subcase  r  =  4  .  The  matrix  is 

0  I-  0 

B4  =  1  0  i 

0^0 


The  eigenvalues  of  this  matrix  are 


The  matrices  Aj^  ,  Aj^^  are  found,  using  standard  techniques,  to  be 

-111 

Aj^  =  0  ^/3 

2  11 


-2  0  2 

^4  =  I  2  n/3  1 

2-^3  1 


The  matrix  E(^(i)  is 

E^Ct) 


4o 


The  sura  of  all  elements  in  is 

The  contribution  to  n(t)  from  subcase  r  a  4  is 


V 


Next  consider  subcase  r  =  3  •  The  matrix  is 


^3  = 


0  i 
1  0 


The  eigenvalues  of  are 


X  =1.  X  a  -  ^ 


The  matrices  ,  A^^  , 


A3. 


V2  -42 


_  k 
*“  2 


1 


1 

42 


_1 

42 


E3(t)  = 


42 


■i/ 


The  sum  of  all  elements  in  A2E3(t)A”^  is 

The  contribution  to  n(t)  from  subcase  r  =»  3  is  therefore 
-J 


noe 


ki 


i 


(2-61) 


(2-62) 


The  last  subcase  Is  r  »  2  .  This  case  is  simply  the  solution  to  (2-36) 


with 

03^ (t)  =  C(t)  . 

The  contribution  is 

nQs""^  .  (2-63) 

Putting  together  the  contributions  [(2-6l),  (2-62)  and  (2-63)]  of  the 
subcases  we  get  the  number  of  molecules  per  unit  volume  with  trapped 
electrons  at  time  t  . 


n(t)  =  n^e”  1  + 


exp  -r  J  +  1 - -  exp  - 


9  +  ^3  /nTs  \ 

-■■■■-  exp  —  J 

6  V  2  / 


exp^ - jj  ,  (2-64) 


J  =  jT C(t')dt'  . 

T 

Recall  that  C(t)  is  the  transition  frequency  for  an  electron  to  escape 
from  a  trap.  Two  expressions  for  c(t)  derived  by  Kubo  are  presented 
as  Eqs.  (A-8)  and  (A-9)  in  Appendix  I.  To  decide  which  is  pertinent  for 


trapping  in  polypeptides,  we  need  to  know  the  vibrational  frequency  of 
the  skeletal  atoms  of  the  molecule.  Infrared  absorption  data^  suggest 

13  -1 

that  the  frequency  is  roughly  2  x  10  sec  .  This  corresponds  to 

_2 

vibrational  quantum  of  hv  9’10  ev  .  Thus  at  liquid  nitrogen  tempera¬ 
ture,  hv/kT  =  10  and  the  first  expression  (A-8)  is  correct .  We  have 

C(t)  =  7  exp^- 


^P.  Doty  and  E.  P.  Geiduschek,  "Optical  Properties  of  the  Proteins" 
in  Proteins,  edited  by  H.  Neurath  and  K.  Bailey,  Vol.  1,  part  A  (Academic 
Press,  Inc.,  New  York,  1953)- 
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The  thermoluminescence  intensity  is  proportional  to  -dn/dt  so  that 


(2-65) 

with  ^ 

J(t)  =  y  Jexp(-  dt')  . 

T 

Equation  (2-65)  represents  the  glow  cvirve  from  a  polymer  with  one  recom¬ 
bination  site  and  three  charge -trapping  sites.  It  may  be  compared  with 
(2-38)  which  represents  the  glow  curve  from  a  molecule  having  one 
recombination  site  and  one  electron-trapping  site.  Such  a  comparison  Is 
made  graphically  in  Fig.  25. 

To  summarize,  In  order  for  a  substance  to  exhibit  thermoluminescence, 
it  or  its  constituent  molecules  must  possess  metastable  excited  states 
which  have  two  properties,  (l)  The  transition  rate  from  the  metastable 
excited  state  to  a  short-lived  excited  state  must  increase  markedly  with 
temperature.  (2)  The  transition  from  this  excited  state  to  the  ground 
state  must  be  accompanied  by  the  creation  of  a  photon.  The  glow  curve 
is  proportional  to  the  time  rate  of  change  of  the  number  of  electrons  or 
molecules  in  metastable  excited  states.  The  problem  of  deriving  a 
formula  for  the  glow  curve  is  essentially  one  of  finding  this  ftmction. 

Glow  curve  derivations  have  been  made  in  terms  of  three  different 
models,  (l)  The  semiconductor  model  involving  valence  and  conduction 
bands  with  traps  for  electrons  and  holes  in  the  energy  gap.  (2)  The 


monomolecular  model  in  which  the  molecules  of  the  solid  have  a  metastable 
excited  state  whose  lifetime  decreases  strongly  as  the  temperature  of  the 
sample  is  increased.  (3)  The  intramolecular  charge-migration  model.  In 
this  model  the  mean  time  for  steps  in  the  migration  process  decreases 
strongly  as  the  temperature  increases.  The  first  derivation  uses  some 
points  of  view  and  results  of  Halperin  and  Braner.^  The  second  derivation 
has  been  made  by  a  number  of  workers.  The  third  is  believed  to  be  a 
new  result. 


■^A.  Halperin  and  A.  A.  Braner,  Phys.  Rev.  117,  ^08  (1960) . 
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CHAPTER  III 


RESULTS  AND  IHTERPRBTATION  OF  THERMOLUMINBSCENCB  EXPERIMENTS 
Thermoluminescence  measurements  have  been  made  in  the  following  way- 
The  sample  -  typically  100  mg  of  powder  or  1  ml  of  liquid  -  is  placed  on 
a  copper  sample  plate  and  covered  with  a  quartz  cover  glass .  A  copper 
constantan  thermocouple  is  attached  to  the  sample  plate  and  the  plate  is 
fastened  to  the  heater.  The  heater  is  mounted  inside  a  Dewar  and  is 
cooled  to  77°  K  by  pouring  liquid  nitrogen  into  the  Dewar  -  see  Figs.  10 
and  11.  Dry  nitrogen  gas  is  fed  into  the  Dewar  to  flush  out  water  vapor, 
thus  preventing  moisture  condensation  on  the  quartz  cover  glass.  After 
the  sample  is  cooled  to  77°  K  it  is  irradiated  with  ultraviolet  light 
from  a  Hanovia  model  306OO  lamp.  There  is  a  band  pass  filter,  consisting 
of  2  Corning  glass  filters,  CS  numbers  a'^d  between  the  ultra¬ 

violet  lamp  and  the  sample.  The  transmission  range  of  this  filter  is 

A  9  XQ 

to  3900  A.  After  the  sample  has  received  10  ergs  of  radiant 

energy,  the  lamp  and  filter  are  removed  and  a  photomultiplier  tube  is 
placed  over  the  sample.  The  heater  current  is  turned  on  and  the  photo¬ 
multiplier  current  and  sample  temperature  are  measured  as  a  function  of 
time. 

Figure  12  is  a  block  diagram  of  the  apparatus  used  for  thermolmines- 
cence  experiments .  Figure  13  is  a  photograph  of  the  cryostat  beneath  the 
photomultiplier  tube  enclosure.  Figure  l4  is  a  photograph  of  the  entire 
apparatus .  The  high  voltage  supply  is  a  Baird  Atomic  power  supply,  model 
3I2A.  It  supplies  voltage  to  a  DuMont  No.  6292  photomultiplier  tube. 

The  light  pipe  is  a  4"  length  of  polished  Incite  rod,  l\"  in  diameter. 

The  Keithley  model  I50  AR  microvolt  ammeter  amplifies  the  current  output 
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FIG.  12--A  block  diagram  of  the  thermoluiainescence  apparatus 

FIG.  13--Side  view  of  the  cryostat  in  position  under  the  black  box 

containing  the  photomultiplier  tube.  One  end  of  the  lucite 
light  pipe  leading  to  the  photomultiplier  tube  may  be  seen 
over  the  cryostat  opening. 
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of  the  photomultiplier  tube  and  drives  the  B22  input  channel  of  the 
Varian  Associates  model  G22  dual  channel  recorder.  The  emf  of  the  sample 
thermocouple  is  balanced  against  the  emf  from  the  crushed  ice  and  water 
reference  junction  and  fed  to  the  A-22-A  input  channel  of  the  chart 
recorder.  This  Input  Is  adjxisted  so  that  the  recorded  thermocouple  emf 
agrees  with  that  measured  with  a  Rubicon  model  2732  Potentiometer.  The 
thermocouple  emf  at  liquid  nitrogen  temperature  agrees  with  the  published 
values^  to  within  0.5^^. 

Voltage  across  the  sample  heater  is  regulated  by  a  Variac .  The 
sample  materials  were  obtained  from  the  California  Corporation  for 
Biochemical  Research  (CCBR)  or  from  Nutritional  Biochemicals  Corporation 
(NBCo)  except  where  otherwise  indicated. 

Data  as  supplied  by  the  chart  recorder  consist  of  the  current  output 
of  the  photomultiplier  tube  and  the  balanced  thermocouple  output  plotted 
as  a  function  of  time.  From  these  data  a  plot  of  the  luminescence  and 
the  temperature  as  a  function  of  time,  can  be  obtained.  Figure  15  shows 
such  a  plot  made  from  the  thermoluminescence  data  for  L-tyrosine.  This 
glow  curve  can  be  analyzed  in  terms  of  the  simple  monomolecular  kinetics 
developed  In  the  previous  chapter.  To  review  briefly,  the  basic  assump¬ 
tion  was  that  each  molecule  of  the  sanple  acted  independently  in  storing 
and  releasing  energy.  It  followed  that  if  the  probability  per  unit  time 
for  emission  of  a  photon  were 

^(t)  =  7  exp[- 


^F.  A.  Mauer,  "Low  Temperature  Equipment  and  Techniques"  in  Format ion 
and  Trapping  of  Free  Radicals,  edited  by  A.  M.  Bass  and  H.  P.  Broida, 
(Academic  Press,  New  York,  i960) . 
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then  the  thermolumineecence  intensity  as  a  function  of  time  would  be 


The  time  t  =  0  may  be  chosen  to  correspond  to  the  beginning  of  the 
glow  curve.  Then  during  the  initial  rise  in  thermoluminescence  intensity, 
the  integral  in  (2-26)  will  be  small,  the  exponential  containing  it  will 
be  essentially  unity,  and  (2-26)  can  be  approximated,  after  taking 
logarithms  as 

£ 

ln[l(t)]  s  constant  -  . 

Figure  16  shows  a  logarithmic  plot  of  the  intensity  of  the  beginning 
of  the  tyrosine  glow  curve  as  a  function  of  the  reciprocal  of  the  tempera¬ 
ture  in  degrees  Kelvin.  The  activation  energy  obtained  fr-'m  the  slope  of 
the  line  drawn  through  these  points  is  0.155  ±  0.015  ev.  The  peak 

temperature,  T  ,  is  124°  Kj  the  heating  rate  at  the  peak  is  26°  K  min”^. 

P 

Putting  these  numbers  into  Eq,  (2-28)  we  get  the  frequency  factor 
y  =  1.3  •  10^  sec"^  . 

A  number  of  amino  acids  and  other  molecules  closely  related  to  them 
have  been  investigated  and  the  thermoluminescence  data,  when  obtained, 
analyzed  as  shown  above.  Table  I  contains  the  averaged  results  of  these 
experiments.  Tyrosine  and  phenylalanine  which  are  discussed  in  detail 
in  this  chapter  are  included  in  the  table  for  the  sake  of  completeness . 
Samples  of  the  same  compound  from  different  biochemical  supply  houses 
show  in  some  cases  different  responses.  The  source  of  the  sample  has 
therefore  been  Included  in  the  table. 

The  glow  curve  of  CCBR  L-phenylalanlne  (Plg.  17)  has  two  peaks  and 
a  strong  initial  luminescence  or  phosphorescence  subsequent  to  its 
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TABLE  1 


This  table  lists  the  amino  acids  and  closely  related  molecules 
which  have  been  investigated  for  thermoluminescence.  All  samples  were 
in  polycrystalline  form.  All  were  cooled  to  liquid  nitrogen  temperature, 
irradiated  with  uv  and  subsequently  heated  and  looked  at  for  evidence 
of  thermoluminescence  as  described  in  the  text.  The  first  column  of 
this  table  gives  the  name  of  the  compound,  the  source,  the  side  group 
or  the  structure  if  the  compound  is  not  an  amino  acid.  The  second 
column  gives  the  maximum  thermoluminescence  intensity  in  units  of 
milllmieroamps,  the  observed  photomultiplier  current.  Column  3  gives 
the  temperature  at  the  peak,  or  peaks  if  there  is  more  than  one  and 
they  are  reasonably  well  resolved.  Columns  4  and  5  give  the  activation 
energy  and  frequency  factor,  determined  in  the  way  described  in  the 
text  whenever  this  has  been  possible. 
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TABLE  X  (Continued) 


COMPOUiro,  SOURCE  AND 
STRUCTURE  OR  SIDE  GROUP 

TL  INTENSITY 
IN  mp  AMP 

T  °K 

P 

E  ev 

y  min  ^ 

L-Leucine  (CCBR) 

CH-CH„- 
/  2 

CH3/ 

<  1 

• 

• 

L-Isoleucine  (CCBR) 
CH^-CHg-CH- 

cii^ 

<  1 

L-Phenylalanine  (CCBR) 

38 

113 

.086+ ,015 

4.10^ 

Q-ch,. 

29 

155 

.221 .06 

.>.5.10’^ 

L-Phenylalanine  (NBCo) 
(Same  side  group) 

60 

123 

.111.015 

8.10^ 

D-Phenylalanine  (CCBR) 
(Same  side  group) 

<  1 

- 

- 

- 

D-Phenylalanine  (NBCo) 

(Same  side  group) 

<  1 

• 

L-Tyroslne  (CCBR) 

ho-Q-ch,. 

3^0 

124 

.1551.015 

8.10^ 

L-Cystelne .HCl  (CCBR) 
HS-CHg- 

<  1 

“ 

L-Methionine  (CCBR) 
CH3-S-CH2-CH2 

<  1 

- 

- 

- 

L-Tryptophan  (CCBR) 

9 

115 

TL  too  weak  for  analysis; 
peaks  were  not  well 

OO'"-- 

1 

H 

10 

152 

separated. 
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TABLE  1  (Continued) 


COMPOUND,  SOURCE  AND  TL  INTBNSITy  „ 
STRUCTURE  OR  SIDE  GROUP  IN  mu  AMP  p 


E  ev 


7  min 


L-Proline  (CCBR)  1  115*5  TL  too  weak  for  analysis 


H_C  CH-COOH 

1 

H 

L-Hydroxyproline  (CCBR)  <  1 

H 

I 

HO-C - CH. 

I  1  " 

H-C  CH-COOH 

i 

L-Histidine  (CCBR)  <  1 

HC=C-CH„- 

I  I  " 

N  NH 

I 

H 

L-Aspartic  Acid  (CCBR)  <  1 

H0-C-CH_- 

II  2 

0 

L-Glutamic  Acid  (CCBR)  <  1 


-v  2  138*5  TL  too  weak  for  analysis 


HO-C-CHg-CHg- 

0 


L-Arginine  (CCBR) 
H2N-J-NH-(CH2)3- 
HH 


L-Lyslne  (CCBR) 
H,N-(CH„),- 


5T 


TiSLS  1  (Continued) 


COMPOUND,  SOURCE  AND  TL  INTENSITY 
STRUCTURE  OR  SIDE  GROUP  IN  mu  AMP 


Phenylglyclne  (NBCo) 


Tyramine-mono  HCl 

(ccbr) 


NHg 


166±5 

194±5 


.08±.015  9.10-^ 

Peak  overlap  prevents 
analys is 

TL  too  weak  for  analysis 
TL  too  weak  for  analysis 


D-Dlhydroxyphenylalanine 

(NBCo) 


TL  too  weak  for  analysis 


3-Hydroxytyranilne  .HCl 
(NBCo) 


122±5  TI>  too  weak  for  analysis 


Indole-3-Acetic  Acid 
(NBCo) 


-CHg-C-OH 

0 


155±5 


TL  too  weak  for  analysis 


This  glow  curve  looks 
like  the  sum  of  several 
unresolved  peaks . 


Tryptamlne.HCl  (NBCo) 


CCBR  -California  Corporation  for  Biochemical  Research 
NBCo  -Nutritional  Biochemicals  Corporation 


irradiation  with  ultraviolet  light  at  77°  K.  A  sanilogarithmic  plot  of 
this  luminescence  shows  that  it  can  be  resolved  into  the  sum  of  two 

exponentied  decays,  one  with  a  half  life  of  3*1  seconds,  the  other  with 

e,  half  life  of  92  seconds.  This  analysis  is  shown  in  Pig.  l8.  Figure  19 
sb-'ws  the  initial  curve  analysis  of  the  first  peak  of  the  L-phenylalanine 
glow  curve.  The  points  seem  to  form  a  tilted  S-shaped  curve  rather  than 

the  straight  line  drawn  through  them.  However  the  deviation  from  a 

straight  line  is  everywhere  less  than  the  error  due  to  the  absolute 
uncertainty  of  the  temperature  (±3°  K)  and  may  be  due  to  systematic  errors 
of  the  temperature  measuring  equipment .  The  error  bars  in  Fig .  19  are 
derived  from  the  absolute  uncertainty  of  the  temperature  measiirements . 

The  overlap  of  the  two  peaks  makes  the  analysis  of  the  second  peak  for 
values  of  E  and  y  quite  xmcertaln. 

Additional  evidence  regarding  the  association  of  monomolecular 
kinetics  with  amino  acid  glow  curves  comes  from  the  comparison  of  the 
theoretical  expression  (2-26)  with  the  experimental  data.  The  Burroughs 
Algebraic  Compiler  at  the  Stanford  Computation  Center  was  used  to  obtain 
numerical  values  for  this  expression.  In  doing  so,  experimental  values 
of  J!  and  y  ,  and  temperature  as  a  function  of  time,  T(t)  ,  were  used. 
The  computer  was  programed  to  evaluate  (2-26)  over  a  time  span  covering 
the  glow  curve,  in  increments  of  .02  minutes.  The  actxial  program  used, 
in  the  language  of  the  computer,  is  given  in  Appendix  II. 

For  the  case  of  L-phenylalanlne,  the  computed  curves  are  compared  to 

the  experimental  data  from  CCBR  L-phenylalanine .  Analysis  of  the  first 

peak  of  a  typical  glow  curve  yielded  the  values  E  =  0.095  ev  , 
h  -1 

7  =  3*ol  *10  min  .  The  theoretical  curves  obtained  using  these  and 
certain  other  values  of  E  ,  7  are  shown  in  Fig.  20 j  the  experimental 
curve  is  also  shown  for  comparison. 
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FIG.  l8--Seinilogarithriic  plot  of  the  phosphorescence  decay 
from  uv  irradiated  L-phenylalanlne  at  77°  K. 

6o 


FIG.  19 — Semllogarithniic  plot  of  intensity  versus  IOO/T  for  the 
initial  part  of  the  phenylalanine  glow  curve. 
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FIG.  20--A  comparison  of  experimental  points  and  theoretical  curves 

computed  using  the  parameters  listed.  The  experimental  points 
shown  as  dots  -  are  from  the  low  temperature  peak  of  an 
L-phenylalanine  glow  curve.  The  solid  lines  are  plots  of 
Eq.  (2-26). 
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ITheoretloal  and  experimental  glow  ciunres  of  L-tyroslne  (CCBR)  are 
shown  in  Pig.  21,  Using  the  twiperature  and  intensity  curves  of  a  typical 
experiment,  we  obtained  the  values  E  -  .158  ev,  7  ■  8.11  •  10^  min"^  . 
These  values  and  the  temperature  curve  were  then  used  In  obtaining  a 
numerical  plot  of  Eq,  (2-26)  with  gn(o)  set  equal  to  1  .  Numerical 

plots  were  also  made  using  the  values  E  =  .I8I  ev,  7  =  8I.6  •  10  min 

^  "2. 

and  E  =  .129  ev,  7  =  .453  •  10  min"  .  The  values  of  7  were  chosen 
so  that  the  curves  peaked  at  the  experimentally  observed  temperature. 

There  Is  a  rather  conspicuous  disagreement  between  experiment  and  theory 
on  the  high  temperature  side  of  the  glow  curve  in  that  the  measured 
thermoluminescence  persists  somewhat  longer  than  theory  would  predict. 
Possible  explanations  are  (l)  the  frequency  factor  7  is  not  a  constant 
with  respect  to  temperature  and  (2)  there  is  a  second  weaker  thermolumines¬ 
cence  peak  which  cannot  be  resolved  but  which  has  a  peak  temperature 
greater  than  125®  K. 

Another  experiment  concerning  the  pertinence  of  the  monomolecular 
model  is  the  following.  A  small  amount  of  L-phenylalanine  was  dissolved 
in  distilled  water  and  one  ml  of  the  solution  frozen  to  77®  K  on  the 
sample  plate.  This  frozen  solution  produced  weak  thermoluminescence 
after  being  irradiated  with  viltraviolet;  a  one  ml  sample  of  frozen 
distilled  water  did  not.  The  thermoluminescence  observed  from  the  frozen 
solution  corresponds  roughly  to  the  low  temperature  peak  of  powdered 
L-phenylalanine.  In  addition,  samples  of  L-phenylalanlne,  D- phenylalanine 
and  L-tyroslne  have  been  freeze  dried  from  distilled  water  solution  and 
found  to  exhibit  thermoluralnescence .  In  fact,  samples  of  D- phenylalanine 
do  not  show  therraolumlnescence  xuitll  freeze  dried  .  These  experiments 
show  that  thermoluminescence  still  occurs  when  the  materials  are  in  a 
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FIG.  21--A  comparison  of  experimental  points  from  an  L-tyrosine 
glow  curve  with  theoretical  glow  curves  computed  using 
the  parameters  listed. 
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dissolved  or  amorphous  state  and  thus  that  crystalline  structure  is  not 
essential. 

Knowledge  of  the  spectrum  of  the  thermoluminescence  would  also  be 
helpful  in  understanding  the  physics  of  the  process.  Unfortunately  the 
intensity  is  too  weak  to  permit  the  luminescence  to  be  passed  through  a 
monochromator  and  still  be  detected.  A  crude  spectrum  has  been  obtained 
in  the  following  way. 

It  was  first  determined  that  repeated  cooling,  ultraviolet  irradiation 
and  warming  of  the  sample  did  not  diminish  the  thermolmlnescence 
response.  Then  a  series  of  experiments  was  performed  with  the  same 
sample,  each  time  covered  with  a  different  Corning  sharp  cut  glass  filter. 
Comparison  of  the  intensity  measured  in  these  experiments  permits  a  rovigh 
estimate  of  percentage  of  luminescence  with  wavelength  greater  than  certain 
values.  The  results  of  these  measurements  are  given  in  Table  2.  The 
filters  used  have  CS  numbers  3-70,  3-72,  3-7^  and  3-75* 


TABLE  2 


Compound 

Percentage  of  Thermoluminescence  with  Wavelength 

(all  CCBR) 

>6000  A 

>5030  A 

>4500  A 

>4080  A 

>3980  A 

L-tyroaine 

0 

16 

37 

73 

95 

L-phenylalanine 
(1st  peak) 

0 

16 

4o 

75 

85 

L-phenylalanine 
(2nd  peak) 

0 

16 

46 

80 

90 

Egg  Albumin 

10 

50 

95 

(Protein) 


The  figures  given  in  this  table  are  estimated  to  be  within  20$  of 
the  correct  values . 
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Dose  response  measurements  have  teen  made  on  L-tyrosine  and  L-phenyl- 
alanine  hy  irradiating  the  sao?>leB  for  varying  lengths  of  time  and  com¬ 
paring  the  relative  6ai?)litudes  of  resulting  thermoluminescence.  For  "both 
confounds,  a  dose  of  2  •  lo"^  ergs  (30  seconds  illumination)  was  suffici- 
ent  to  produce  the  maximum  level  of  response,  while  a  dose  of  6  •  10° 
ergs  produced  half  the  maximxim  level.  The  quantum  yield  of  the  thermo- 

O 

luminescence  produced  hy  these  materials  is  roughly  one  part  in  10°. 

Next  to  he  considered  is  the  nature  of  the  metastahle  excited  states. 
The  three  hypotheses  proposed  in  Chapter  II  were  (l)  trapping  at  an  ionic 
portion  of  a  molecule,  (2)  a  frozen  isomerization,  and  (3)  a  photo  induced 
charge -transfer  conplex  involving  a  metal  ion.  The  evidence  which  is 
helpful  in  deciding  which  of  these  is  correct  is  (a)  only  about  one 

Q 

molecule  in  10  develops  this  metastable  state,  (b)  the  same  compound 
from  different  biochemical  suppliers  shows  different  thermoluminescence 
responses,  (c)  the  naturally  occurring  L  isomer  of  phenylalanine  develops 
thermoluminescence,  while  the  artificially  synthesized  D  isomer  does  not. 
These  three  pieces  of  evidence  seem  most  consistent  with  the  third 
hypothesis.  It  appears  that  trace  amounts  of  some  metal  ion  come  through 
the  Isolation  procedures  with  the  naturally  occurring  amino  acids  and 
proteins,  but  do  not  occur  in  the  laboratory  synthesized  D  isomer  of 
phenylalanine.  The  actual  presence  of  trace  amounts  of  metal  ions  has 
been  determined  in  two  ways.  Distilled  water  solutions  of  L-tyrosine  and 
L-phenylalanine  were  shaken  with  a  solution  of  the  chelating  agent 
8-qulnollnol  in  CCl^^.  After  a  few  hours,  a  distinct  color  appeared 
indicating  the  presence  of  copper  or  iron  ions  or  both.*  Wide  sweep 


*I  am  indebted  to  L.  Taskovich  for  this  technique. 
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BPR  measurement B  of  tyrosine  and  phenyleileHiine  samples  revealed  the 
presence  of  paramagnetic  Ions,,  their  abundance  being  a  few  parts  per 
million.*  Further,  tyrosine  was  dissolved,  filtered  and  recrystallizied 
from  water  solution.  .The  thermoluminesoence  response  of  the  purified 
material  was  about  of  the  untreated  tyrosine.  The  filtrate  was 
freeze  dried  and  found  also  to  exhibit  thermoluminescence.  It  appears 
that  some  substance  which  was  present  in  the  untreated  sample  was 
partially  removed  by  the  recrystalllzatlon,  going  at  least  in  part  into 
the  filtrate.  This  result  suggested  the  following  experiment.  It  seemed 
likely  that  if  an  ion,  present  in  trace  amounts,  were  necessary  for 
thermoluminescence,  adding  an  abundance  of  this  ion  would  greatly  increase 
the  intensity  of  the  response.  Table  3  summarizes  the  results  of  adding 
one  part  of  the  salts  listed  to  a  water  solution  of  ten  parts  of  (CCBR) 
L-phenylalanine,  freeze  drying  the  mixture  and  examining  for  thermo¬ 
luminescence.  Of  the  ten  salts  tried,  only  two  enhance  the  response, 
only  one  appreciably.  Copper  and  iron  salts  quench  the  response  entirely. 
Thus  copper  and  iron  are  removed  from  suspicion  as  the  unknown  coraplexed 
quantity  and  calcium  is  introduced  as  the  only  likely  agent  of  the 
materials  tested.  Further  work  and  more  sophisticated  chemistry  are 
needed  before  the  amino  acid-metal  ion  complex  hypothesis  can  be  accepted 
with  certainty. 

Polymers  of  amino  acids  have  been  found  to  exhibit  thermoluminescence 
althouigh  the  results  are  not  what  one  would  predict  from  the  monomer  case. 
Poly  L-tyrosine,  for  example,  shows  a  weak  response  and  poly  L-phenylalanine 
a  stronger  response  but  only  one  peak.  (See  Fig.  22.)  These  results 


*The8e  measurements  were  kindly  performed  by  Dr.  J.  Mallng. 


TABLE  3 


Salt  Added 

Result 

BaCl 

Intensity  of  first  peak  normal 

Intensity  of  second  peak  reduced  to  half  value 

MgSOi^ 

Intensity  of  first  peak  normal 

Intensity  of  second  peak  reduced  to  about  half 
value 

MnClg 

Intensity  of  both  peaks  reduced  to  about  one-third 
normal  value 

CoClg 

Intensity  reduced  to  about  one-fourth  normal  value 

WagMoOj^-aH^O 

Intensity  of  both  peaks  reduced  to  about  one-third 
normal  value 

ZnSOj^ 

Intensity  of  first  peak  enhanced  by  about 

Intensity  of  second  peak  about  half  normal  value 

FeClg 

This  salt  completely  quenches  both  peaks 

CaClg 

Intensity  of  first  peak  increased  by  a  factor  of  3 
Second  peak  about  one-foiurth  normal  value 

CUClg-SHgO 

Both  peaks  completely  quenched 

CuCl 

Both  peaks  completely  quenched 
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and  the  results  of  experiments  on  di-  and  tripeptides  and  on  proteins 
are  summarized  in  Table  4. 

The  iron-free  proteins  which  were  examined  showed  a  very  long-lived 
luminescence  at  8o°  K.  The  decay  is  complex  hut  can,  in  the  case  of  egg 
albumin,  be  resolved  as  shown  in  Fig,  sU,  The  lifetime  of  this  lumines¬ 
cence  is  very  temperature  dependent.  Had  lower  temperatures  been  avail¬ 
able,  this  luminescence  could  have  been  frozen  in,  with  much  longer 
lifetimes.  Because  the  luminescence  leaks  out  at  8o°  K,  the  intensity 
of  the  peak  varies  considerably.  The  glow  curve  of  trypsin  was  found  to 
be  more  reproducible j  it  is  shown  in  Fig.  23. 

Note  that  this  curve,  in  contrast  to  those  of  Figs.  15  and  17  has  a 
steeper  slope  on  the  ascending  than  the  descending  side.  This  contrast 
is  consistent  with  a  change  from  monomolecular  to  polymeric  kinetics. 
Figure  25  shows  graphs  made  of  (2-26)  and  (2-65)  using  the  same  values 
of  E,  y  and  the  same  dependence  of  T  upon  t.  Both  curves  were  normal¬ 
ized  by  setting  the  integral  of  the  intensity  equal  to  1.  The  monomer 
curve  has  its  steepost  slope  on  the  descending  side,  while  the  L  =  4 
polymer  curve  is  nearly  symmetric  but  is  slightly  less  steep  on  the 
descending  side.  This  asymmetry  will  be  pronounced  for  larger  L. 

Thus  we  conclude  that  the  polymer  glow  curves  are  consistent  with 
the  charge  migration  hypothesis  and  that  the  process  is  an  intramolecular 
one.  The  absence  of  thermolumlnescence  from  ferritin  and  cytochrome  C 
and  also  from  a  lyophylized  mixture  of  egg  albumin  and  ferris  chloride 
indicates  that  iron  somehow  quenches  the  luminescence  from  polypeptides 
as  well  as  from  L-phenylalanlne . 
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TABLE  L 


Compound 


Olycyl  L-Phenylalanine 
NBCo 

Glycyl  L-Tyroaine 
HBCo 


DL  Leucyl -Glycyl 
DL  Phenylalanine 
NBCo 

Poly  L-Tyroslne 
Yeda 

Poly  L-Phenylalanine 

Trypsin 

CCBR 

Egg  Albumin:  10  mg  in 
1  ml  distilled  H^O 
CCBR 


Chymotrypsin 

CCBR 

Cytochrome  C 
CCBR 

Ferrlt in 
CCBR 


Results 


One  peak,  intensity  8  mp  amp,  =  108°  K 

One  peak,  intensity  12  mp  amp,  Tp  =  136°  K, 
a  shoulder  on  low  temperature  at 
T  =  104°  K 

One  peak,  intensity  7  md  amp,  Tp  =  116°  K 
with  shoulders  on  both  high  and  low 
temperature  sides 

One  peak,  intensity  10  rafi  amp,  long  tall 
on  high  temperature  side 

See  Pig.  22 

See  Fig.  23 


One  peak,  intensity  about  50  mp  amp, 
but  dependent  upon  length  of  time 

after  uv  irradiation  -  T  =  120°  K 

P 

Sample  shows  long-lived  decay  and  glow 
curve  much  like  that  of  Trypsin 

No  thermoluminescence  or  long-lived 
initial  glow 

No  thermoluminescence  or  long-lived 
Initial  glow 


THERMOLUMINESCENCE  OF  TRYPSIN 
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MINUTES 


INTENSITY  (m /I  amps) 


FIG.  24 — Semllogarithmlc  plot  of  the  decay  of  phosphorescence 
from  uv  irradiated  egg  albumin  at  77°  K. 


RELATIVE  INTENSITY 


FIG.  25--Comparlson  of  theoretical  glow  curves  from  a  monomer 
and  a  four-unit  polymer.  The  same  values  of  E  ,  y 
and  T(t)  were  used  in  computing  these  curves. 
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CHAPTER  IV 


IQBULT3  AHD  tM'Ma>m!ATION  OF  m9CTR0H 
FARAMAOMriC  B^NANOE  BXFimiMENTS 

When  polycryatalllne  samples  of  amino  acids,  polypeptides  or  proteins 
are  irradiated  with  ultraviolet  light  at  80®  K,  they  develop  paramagnetism 
which  can  be  measured  by  means  of  an  EER  (electron  paramagnetic  resonance) 
spectrometer.  This  paramagnetism  can  be  interpreted  as  arising  from  free 
radicals  or  trapped  electrons  or  ions  -  products  of  the  absorption  of  the 
ultraviolet  light  by  the  molecules  of  the  sample. 

The  physics  of  EER  absorption  by  unpaired  spins  is  fully  discussed 
in  the  literature.^  For  the  studies  reported  here,  one  may  describe  the 
phenomenon  In  terms  of  a  dilute  ensemble  of  unpaired  electrons  in  the 
material  under  invest Igatipn.  These  electrons  interact  (l)  with  an 
externally  applied  magnetic  field,  (2)  with  a  rapidly  varying  magnetic  field 
in  the  microwave  cavity  containing  the  material,  (3)  with  the  electro¬ 
static  crystal  field  due  to  surrounding  atoms  or  molecules,  {k)  with  the 
magnetic  moments  of  neighboring  atomic  nuclei. 

The  first  two  interactions  may  be  represented  by  the  spin  Hamiltonian 
X  =  g^HQS^  +  gp2H^co8(urt)S^  (J^-1) 

where  g  is  the  spectroscopic  splitting  factor 
P  is  the  Bohr  magneton 

Hq  is  the  externally  applied  magnetic  field  which  defines 
the  z-axls 

is  the  z-component  of  the  spin  operator 

^G.  E.  Eake,  Earamagnetlc  Resonance  (Benjamin,  New  York,  I962) . 


2H^co8(ut)  Is  the  microwave  magnetic  field 
is  the  x-component  of  the  spin  operator. 

The  second  term  in  (4-l)  serves  as  a  perturbation  which  causes  transitions 
between  eigenstates,  ■^  >  =  ,  of  the  first  term.  The  probability  of 

these  transitions  is  proportional  to 


8ln®(V2)t 
.a  )  A® 


(4-2) 


where 


tt 


There  are  two  effects  which  cause  the  smearing  out  of  the  resonance 
condition, 

=  hu  .  (4-3) 

(1)  Because  of  the  electron-crystal  field  interaction,  the  g-factor  is 

slightly  dependent  upon  the  orientation  of  the  crystal  field  relative 
to  Hq  .  In  polycrystalline  samples  where  all  orientations  occur,  g 
has  a  range  of  values,  with  <!^/g  10  in  amino  acid  samples. 

(2)  There  are  local  variations  in  within  the  sample.  These  two 
effects  may  be  taken  into  account  by  introducing  a  density  function 

f(u)  for  the  resonance  frequency.  Then  one  can  show  that  the  transition 


rate  between  eigenstates  of  is 


(4-4) 


Experimentally,  the  resonance  is  detected  when  the  transitions  result 
in  the  sample  absorbing  power  from  the  microwave  field.  If  N  is  the 
number  of  unpaired  spins  in  the  sample,  and  W  is  the  spontaneous 
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spin-flip  probability,  the  power  absorbed  Is 

hwVN  tanh(gPHV2l5l) 

P(«) - ^ -  ,  {h-3) 

1  +  2T^V 

where 

1 

T  =  -  . 

W[1  +  exp(gPHQ/kT)] 

The  interaction  between  the  spin  of  the  unpaired  electron  and  the 
moments  of  neighboring  atomic  nuclei  is  called  the  hyperflne  interaction. 

It  can  be  treated  by  adding  the  term 

Nuclei 

to  (4-1 ) .  In  the  representation  in  which  the  first  term  of  (4-1)  is 
diagonal,  (4-6)  combines  with  this  term  to  give 

<H>  s  ggH^m  +  A^M^m  .  (4-7) 

Nuclei 

Here,  is  the  magnetic  quantvun  number  of  the  1th  nucleus,  m  that 

of  the  electron  and  is  the  strength  of  the  interaction.  Equation  (4-7) 
may  be  rewritten  as 

<^>=g^(Ho+^  (4-8) 

showing  that  this  interaction  produces  a  shift  in  the  resonance  condition. 
Consider  this  special  case  where  an  electron  interacts  equally  with  two 
hydrogen  nuclei.  ,  The  sura  in  (4-8)  can  take  on  three  values  depending 
upon  the  orientation  of  the  two  protons'  magnetic  moments  relative  to 
.  The  three  possible  orientations  are:  (l)  both  moments  parallel 
to  ;  (2)  both  moments  antiparallel  to  j  (3)  one  moment  parallel 
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and  one  moment  antiparallel.  Thus,  electrons  interacting  equally  with 
two  protons  will  undergo  spin-flip  transitions  when  u)  »  (case  3) 

and  when  cj  ■  i3[HQ  ±  (2A/gP)]  (cases  1  and  2).  The  resulting  three-line 
absorption  pattern  is  believed  to  accomt  for  part  of  the  observed 
absorption  by  ultraviolet  irradiated  glycine.  (See  Fig.  26.) 

Our  EPR  studies  have  been  directed  toward  the  following  questions. 

(1)  Which  amino  acids,  related  molecules  and  proteins  develop 
paramagnetism  when  irradiated  with  ult-ravlolet  light  at  77°  K? 

(2)  What  is  the  location  of  the  unpaired  electron  which  is  presumably 
responsible  for  the  resonance? 

(3)  Are  unpaired  electrons  producing  the  resonance  also  involved 
in  the  production  of  theimolumlnescence? 

The  following  apparatus  and  materials  were  used  in  the  experimental 
work:  (l)  A  Varlan  Associates  EPR  spectrometer  and  six-inch  electromagnet 

form  the  basic  instrument.  (2)  A  quartz  liquid  nitrogen  Dewar  designed 
to  fit  inside  the  microwave  cavity  is  used  to  keep  the  sample  at  77°  K 
during  ultraviolet  irradiation  and  EPR  absorption  measurements.  Quartz 
tubes  3  nro  or  4  mm  in  diameter  connaln  the  sample  materials.  (3)  A 
Hanovla  type  30600  l40  watt  ultraviolet  lamp  is  used  for  irradiation. 

The  light  passes  through  a  Corning  Vycor  glass  filter  before  reaching 
the  Dewar.  This  filter  absorbs  light  of  wavelength  shorter  than  23OO  A. 

(4)  Biochemicals  used  are  obtained  from  CCBR  and  NBCo  and  are  the  purest 
available  from  those  companies . 

The  experimental  procedure  has  been  to  seal  the  compounds  as  they 
are  supplied  by  the  vendor  inside  clean  quartz  tubes,  to  cool  the  materials 
by  immersing  the  tubes  in  liquid  nitrogen,  to  irradiate  the  tubes  while 
they  are  immersed  in  liquid  Kg  ,  to  measure  the  EPR  absorption  of  the 
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6LYCYLGLYCINE 


FIG,  26— First  derivative  traces  of  paramagnetic  resonance  absorption 

in  glycylglycine  and  glycine  measured  as  a  function  of  magnetic 
field  H  .  The  amplitude  of  the  100  KC  modulation  of  Hq  was 
1.5  gauss  peak  to  peak  with  glycylglycine,  2.5  gauss  with  glycine. 
Maximum  value  of  microwave  magnetic  field  at  sample  was 
Hi  =  0.1  gauss  in  each  case.  These  compounds  were  irradiated 
in  the  dry  state  with  ultraviolet  light  and  the  resonance 
absorption  was  measured  afterward  without  changing  the  sample 
temperature  from  77°  K.  Position  of  the  DPPH  absorption 
(g  =  2.0036)  is  indicated  by  an  arrow. 
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sample  at  77°  K  and  again  at  room  tmperature.  The  apparatus  produces 
a  chart  recording  of  the  first  derivative  of  the  absorption. 

A  nvanber  of  amino  acids  develop  only  a  weak  absorption  when  treated 
in  this  way.  Glycine  is  a  typical  example  of  this  class  of  amino  acids. 
Its  first  derivative  trace  is  shown  in  Fig.  26.  Although  this  signal  is 
weak,  as  indicated  by  the  relative  noise  level,  the  shoulders  indicating 
a  three-line  absorption  are  real  and  reproducible  as  has  been  confirmed 
by  Slmmonds.®  This  three-line  component  of  the  signal  may  be  due  to  a 
free  radical  Intermediate  of  photo-induced  deamination.  This  reaction 
is  common  to  moat  amino  acids^^'^  and  would  be  responsible  for  part  of 
the  EPR  signal  developed  by  them.  In  the  case  of  glycine  this  reaction 
may  proceed  in  part  as  shown  in  (4-9) . 


^COOH 

NHg 


OOH 


.CHg  +  MHg 


(4-9) 


where  the  dot  at  the  a-carbon  represents  an  unpaired  electron. 

The  signal  produced  in  glycylglyc ine  which  shows  no  triplet  character 
indicates  that  some  other  free  radical  must  also  be  produced,  Without 
doing  studies  on  single  crystals,  as  Gordy®  has  done  using  x-rays,  it 
is  not  possible  to  form  an  opinion  about  the  structure  of  this  radical. 
Other  amino  acids  which  develop  this  weak  single  absorption  line  are 
L-glutamic  acid,  L-cystlne,  L-histidlne,  L-tryptophan  and  L-phenylalanlne . 


®J.  A.  Simmonds,  private  communication. 

^J.  Greenstein  and  M.  Winity,  Chemistry  of  the  Amino  Acids  (John 
Wiley  and  Sons,  New  York,  i960). 

^A.  Horsfeld,  J.  Morton,  D.  Whlffin,  Trans.  Faraday  Soc.  1657 

(1961). 

®M.  Katayama  and  W.  Qorc^,  J.  Chem.  Phys .  117  (1961) . 

80 


In  each  case,  warming  the  sample  to  room  tanperature  reduced  the  ETR 
signal  intensity  by  roughly  a  factor  of  five. 

Cysteine,  when  treated  as  described,  developed  a  strong  five-line 
absorption.  This  is  shown  in  Tig.  27.  The  shape  of  this  trace  Is 
consistent  with  the  hypothesis  that  the  unpaired  electron  interacts  about 
equally  with  four  protons,  suggesting  that  irradiation  induces  dimerization 
with  an  unpaired  electron  at  the  disulfide  bond,  as  Indicated  in  the 
sketch  below. 

HOO^  COOH 

CH  -  CH-  -  S  1  S  -  CH„  -  CH 
/  2  2  \ 

HgN 

It  is  further  observed  that  when  the  sample  is  warmed  to  room  temperature, 
the  absorption  changes  from  a  five-line  to  a  three-line  pattern  with  about 
20^6  of  the  low  temperature  intensity.  This  may  be  the  result  of  the 
breaking  apart  of  the  free  radical  dimer  and  the  formation  of  a  cysteine 
free  radical  with  the  unpaired  electron  at  the  S  -  CH^  region.  Cystine, 
the  disulfide  bond  dimer  of  cysteine  develops  only  the  weak  single-line 
absorption  shown  at  the  top  of  Fig.  27,  when  treated  as  described.  EPR 
studies  of  cystine  following  x-ray  irradiation®  have  shown  that  the  sulfur 
atom  is  able  to  accept  an  unpaired  electron  even  when  this  atom  is 
involved  in  a  disulfide  bond.  The  resonance  observed  looks  much  like  the 
room-tempe nature  trace  of  cysteine  resonance  in  Fig.  27.  The  dimer  cystine 
thus  seems  rather  impervious  to  ultraviolet  light,  though  not  to  x-rays, 
in  marked  contrast  to  the  monomer  cysteine. 


®J.  Henriksen  in  Free  Radicals  in  Biological  Systems  (edited  by 
M.  S.  Blois,  ^al..  Academic  Press,  New  York,  19^1;. 
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L  -  CYSTINE 


V 


FIG.  2T--First  derivative  traces  of  paramagnetic  absorption  in  L-cystine 
and  L-cystelne.  Measurements  on  L-cystlne  were  made  using  2.5 
gauss  peak-to-peak  100  KC  modulation  and  an  amplitude  of 
0.15  gauss.  For  cysteine,  the  100  KC  modulation  amplitude  was 
0.15  gauss  peak-to-peak  at  both  temperatures.  amplitudes 
used  were  0.07  gauss  at  77°  K,  0.10  gauss  at  room  temperature. 
In  each  case  the  compounds  were  ultraviolet  irradiated  in 
polycrystalline  form  at  77°  K. 


Although  glutamic  acid  and  glycine  do  not  develop  a  strong  signal, 
their  presence  with  cysteine  in  the  tripeptide  glutathione  causes  the 
EPR  absorption  developed  by  reduced  glutathione  to  be  much  different  from 
the  sum  of  the  spectra  of  the  constituents.  See  ?ig.  28.  This  absorp¬ 
tion  seems  to  be  the  superposition  of  a  strong  single  absorption  line 
upon  the  cysteine  absorption.  When  a  mechanical  mixture  of  cysteine, 
glycine  and  glutamic  acid  is  irradiated  with  ultraviolet  light,  only  the 
cysteine  absorption  is  observed.  It  thus  appears  that  the  formation  of 
peptide  bonds  makes  available  new  sites  for  unpaired  electron  trapping 
or  free  radical  formation. 

The  EPR  absorption  developed  in  a  number  of  aromatic  amino  acids 
and  related  molecules  has  been  measured.  Two  classes  of  results  are 
observed.  Those  compounds  which  contain  hydroxy  groups  on  the  ring 
develop  a  strong  single  absorption  line.  Examples  are  tyrosine,  tyramlne, 
dlhydroxy-phenylalanlne,  and  3-hydroxy  tyromlne.  On  the  other  hand, 
those  compounds  such  as  L-phenylalanine,  which  lack  an  OH  group  on  the 
ring  show  a  weaker  signal.  The  role  played  by  the  hydroxyl  group  is 
not  clear.  One  possibility  is  that  the  OH  group  -  being  somewhat  acidic 
when  attached  to  a  ring  structure  -  is  dissociated  by  the  ultraviolet 
light  giving  off  the  proton,  leaving  an  unpaired  electron  on  the  oxygen 
after  the  resulting  rearrangement.  A  search  was  made  for  atomic  hydrogen 
absorption  in  a  sample  of  tyramlne.  It  was  not  found.  It  may  be  that 
it  could  escape  even  at  77*^  K  or  that  it  was  not  formed  by  the  ultraviolet 
irradiation. 

Proteins  develop  EPR  absorption  signals  which  do  not  resemble  those 
of  any  individual  amino  acid.  Trypsin  absorption  is  shown  in  Pig.  29. 

Egg  albumin  absorption  at  77°  K  and  at  300°  K  is  shown  in  Fig.  30. 
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FIG.  28— First  derivative  traces  of  paramagnetic  resonance  absorption 

in  reduced  glutathione.  This  compound  is  a  tripeptide  composed 
of  glutamic  acid,  glycine  and  cysteine.  This  compound,  in  the 
form  of  a  dry  powder,  was  placed  in  a  tube  and  the  tube  was 
evacuated  in  an  attempt  to  remove  oxygen.  The  tube  was  then 
sealed  and  the  compound  was  cooled  to  77°  K  and  irradiated  with 
ultraviolet  light.  During  the  absorption  measurements  at  77°  K, 
the  100  KC  modulation  amplitude  was  0.6  gauss  peak-to-peak  and 
was  approximately  0,15  gauss  at  room  temperature,  the  100  KC 
modulation  was  2.0  gauss  peak-to-peak,  the  field  strength 
was  unchanged. 
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TRYPSIN  AT  7  7*K 


TRYPSIN  AT  ROOM  TEMPERATURE 


I40g 


FIG.  29--First  derivative  trace  of  paramagnetic  resonance  absorption 

developed  in  polycrystalline  trypsin  after  viltraviolet  irradi¬ 
ation  at  77°  K.  At  77°  K  the  H3_  field  was  0.07  gauss,  at 
room  temperature  it  was  0.1  gauss. 
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EGG  ALBUMIN  AT  77®K 


EGG  ALBUMIN  AT  ROOM  TEMPERATURE 


PIG,  3n--Plrst  derivative  traces  of  paramagnetic  resonance  absorption 

developed  in  egg  albumin  after  ultraviolet  irradiation  at  77°  K. 
The  H,  field  was  0.1  gauss  at  both  77°  K  and  room  temperature. 


Allen  and  Ingram'^  find  that  egg  albui^-h  develops  an  EPR  absorption  when 
irradiated  with  3^50  A  at  77°  K  which  disappears  completely  when  the 
sample  is  warmed  to  room  temperature.  When  they  irradiated  with  2537  A 
ultraviolet  light  they  produced  an  EJR  signal  which  did  not  decrease  when 
the  sample  was  warmed  to  room  temperature.  Since  the  light  used  in  our 
experiments  contained  both  of  the  wavelengths  used  by  these  workers  it 
may  be  that  the  signal  shown  at  the  top  of  ?ig.  30  is  the  sum  of  the 
3650  A  and  2537  A  response  while  the  bottom  trace  Is  that  of  the  2537  A 
response  alone.  Allen  and  Ingram  interpret  their  results  by  assuming 
that  molecules  of  egg  albumin,  in  the  presence  of  loosely  bound  water 
are  able  to  act  as  band  model  semiconductors .  They  say  "the  most  straight¬ 
forward  interpretation  of  the  results  on  the  egg  albumin  is  that  the 
electrons  are  being  excited  into  the  upper  conduction  band  by  the  3^50  A 
irradiation  and  some  of  these  are  then  falling  into  low-lying  traps 
formed  by  the  presence  of  odd  protons  along  the  polypeptide  chain  system. 
On  warming  to  room  temperature,  these  trapped  electrons  are  re-excited 
to  the  conduction  band  and  then  return  to  the  ground  state  so  that  the 
resonance  signal  disappears."  The  results  obtained  by  these  authors 
could  as  well  be  Interpreted  in  terms  of  the  Markov  process.  The 
experimental  evidence  does  not  indicate  which  of  these  two  hypothesis 
is  closer  to  the  truth.  In  fact  it  is  not  yet  clear  even  in  the  case 
of  the  relatively  simple  anthracene  molecular  crystal  whether  charge- 
carrier  motion  should  be  described  by  the  band  model  or  the  hopping 
model . 


’’^B.  T.  Allen  and  D.J.E,  Ingram,  Compte-rendu  du  9^  Colloque  Ampere, 
p.  219,  Pise,  September  i960. 
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Bovine  hemoglobin  when  irradiated  at  77°  K  develops  the  BHl  absorption 
shown  in  Fig.  31.  The  structure  shown  in  this  trace  has  not  been  inter¬ 
preted,  This  absorption  disappears  below  the  weak  pre irradiation  signal 
when  the  sample  is  warmed  to  room  temperature.  It  is  remarkable  that  the 
presence  of  iron  causes  the  quenching  of  EER  signals  at  room  temperatures 
and  also  of  thermoluminescence  from  proteins. 

Our  EER  studies  have  shown  that  (l)  ultraviolet  Irradiation  at  77°  K 
causes  some  absorption  In  many  amino  acids  and  probably  all  proteins. 

The  ultraviolet  Induced  deamination  reaction  (^-9)  which  results  in  an 
unpaired  electron  at  the  os-carbon  may  be  responsible  for  at  least  part 
of  the  absorption  observed  In  glycine,  L-glutamic  acid,  L-cystine, 
L-histidlne,  L-tryptophan  and  L-phenylalanlne .  L-cysteine  develops  an 
absorption  pattern  which  suggests  that  ultraviolet  irradiation  dimerizes 
this  molecule  at  77°  K.  Aromatic  amino  acids  with  OH  groups  attached 
to  the  ring  develop  a  strong  slngle-llne  absorption  suggesting  the  OH 
group  loses  Its  proton  when  the  molecule  absorbs  an  ultraviolet  photon. 

The  trlpeptlde  glutathione  develops  an  absorption  pattern  which  Is  more 
than  the  sum  of  the  absorption  patterns  of  the  three  constituents, 
glycine,  L-glutamic  acid  and  L-cystelne.  All  of  the  proteins  examined, 
trypsin,  egg  albumin  and  bovine  hemoglobin  develop  absorption  after 
ultraviolet  Irradiation.  There  Is  some  poorly  resolved  hyporflne 
structure  which  does  not  reveal  the  location  of  the  unpaired  electrons 
in  the  sample,  (in  polycrystalline  material  the  smearing  of  the  resonance 
condition  (4-3)  results  in  the  smearing  of  the  hyperflne  structure. 
Conclusions  based  on  interpretation  of  hyperfine  structure  must  therefore 
be  considered  tentative.)  (2)  There  is  no  evidence  that  the  unpaired 
electrons  responsible  for  EER  absorption  are  also  Involved  in  the 


88 


BOVINE  HEMOGLOBIN  AT 


77*K 


FIG.  31 — First  derivative  trace  of  paramagnetic  resonance  absorption 
developed  in  bovine  hemoglobin  after  ultraviolet  irradiation 
at  77°  K.  Only  the  weak  pre- irradiation  absorption  could  be 
detected  in  this  sample  when  it  was  warmed  to  room  temperature, 
field  was  approximately  O.O7  gauss. 
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production  of  thermolumlnescenee .  The  number  of  unpaired  spins  In 
ultraviolet  Irradiated  amino  acids  Is  several  orders  of  magnitude  greater 
than  the  number  of  photons  emitted  during  thermolumlnescenee. 
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SUMMARY 

The  research  reported  in  this  thesis  has  shown  that  certain  amino 
acids,  other  closely  related  molecules,  amino  acid  polymers  and  proteins 
develop  metastahle  excited  states  with  very  long  lifetimes  (hours)  when 
irradiated  with  ultraviolet  light  at  77°  K.  Transitions  from  these  states 
to  states  of  lower  energy  result  in  the  emission  of  light.  This  light, 
in  the  form  of  thermoluminescence,  has  been  observed  upon  heating  the 
compounds  following  ultraviolet  irradiation.  The  wavelength  of  exciting 
light  is  in  the  range  2t-00  S  X.  S  3900  A.  The  maximum  of  the  observed 
themnolumlnescence  is  in  the  region  from  UOOO  to  A5OO  A.  Experimental 
results  are  consistent  with  the  hypothesis  that  the  metastable  excited 
state  is  molecular  rather  than  crystalline  in  nature,  (l)  Thermolumlnes- 
cence  has  been  detected  from  compounds  in  aqueous  solution  and  after 
freeze  drying.  (2)  The  luminescence  intensity  can  be  described  in  terms 
of  monomolecular  kinetics.  (3)  The  wavelength  of  the  absorbing  light  is 
in  the  region  of  an  absorption  band  of  the  molecules  Involved.  The 
metastable  state  Involved  is  believed  to  arise  in  a  charge  trapping 
complex  of  the  molecule  and  a  metallic  ion.  A  new  expression  for  th^rrao- 
luminescence  Intensity  from  polymers,  based  on  a  Markov  process  charge 
migration  model,  has  been  derived. 

We  have  shown  that  ultraviolet  irradiation  of  these  compounds  at  77°  K 
also  produces  in  many  cases  electron  paramagnetic  resonance  absorption. 

The  aromatic  amino  acids  develop  a  single  absorption  line  with  a  g-value 
close  to  that  of  DPPH.  The  intensity  of  absorption  is  greater  in  those 
molecules  which  contain  hydroxy  groups  on  the  ring  than  in  those  which 
^  do  not.  Cysteine  has  been  found  to  develop  a  hyperfine  structure  which 
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suggests  that  a  dimer  is  formed.  The  absorption  spectra  of  several 
other  amino  acids  and  several  proteins  have  been  measured.  In  general 
the  intensity  of  the  absorption  is  much  reduced  after  the  samples  are 
warmed  to  room  temperature.  It  has  not  been  possible  however  to  show  a 
correlation  between  the  EPR  absorption  and  therraolumlnescence .  The 
absorption  is  believed  due  to  trapped  free  radical  Intermediates  of 
photo-chemical  reactions. 

This  research,  particularly  the  thermolumlnescence  work,  supports 
the  belief  that  charge  migration  via  some  hopping  mechanism  may  take 
place  within  protein  molecules.  This  work  shows  that  ultraviolet  light 
may  initiate  the  charge  migration,  and  It  suggests  that  metal  ions  play 
an  Important  role  in  the  charge  migration  and  trapping  processes. 
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Asmoax  I 


The  Born-Oppenhelmer  approximation  can  be  used  to  derive  what  has 
been  called  C(t)  ,  the  probability  per  unit  time  that  a  charge  will  be 
thermally  excited  out  of  a  trap.  The  trap  is  thought  of  as  a  positive 
atom  or  ion  situated  in  a  solid.  The  Hamiltonian  for  the  electron- ion 
system  can  be  written  as 

+  V  (A-1) 

where  ia  the  kinetic  and  potential  energy  operator  for  the  positive 

atom,  represents  the  kinetic  energy  of  the  electron  and  its  inter¬ 

action  with  the  rest  of  the  solid  or  molecule,  and  V  is  the  interaction 
of  the  electron  with  the  positive  atom.  In  the  calculation,  V  rather 
than  the  Interaction  with  the  molecule  or  solid  as  a  whole  is  assumed  to 
determine  the  shape  of  the  trapped  electron-wave  function. 

The  essence  of  the  first-order  Born-Oppenhelmer  approximation  is  the 
factorization  of  the  elgenfvinction  of  (A-l)  into  two  functions  in  the 
following  way.  Considering  the  coordinates  Q  of  the  trapped  atom  and 
the  rest  of  the  solid  as  fixed,  one  finds  the  solution  of  the  Schrodlnger 
equation  for  the  electron  with  coordinates  q  to  be 

[Jig  +  V(q,Q)]<p^(q,Q)  =  €^(Q)cPi(q,Q)  .  (A-2) 

The  eigenvalue  e^(Q)  which  appears  in  this  equation  serves  as  the 
potential  produced  by  the  electron  at  the  position  of  the  positive  atom 
in  the  Schrodlnger  equation  for  this  atom,  which  is 

The  factorization  mentioned  above  consists  of  writing  the  elgenfimction 
of  (A-1)  as 

=  CPi(q,Q)X^g(Q)  .  (A-k) 
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Operating  on  (A-4)  with  (A-l),  one  obtains 

where 

The  term  is  nonvanishing  because  the  differential  operators  in 

Jljj  act  upon  the  atomic  coordinatee  Q  which  appear  as  parameters  in 

cpi(<l,Q)  . 

The  operator  Jl'  is  the  perturbation  which  provides  the  electron- 
phonon  interaction  necessary  for  the  thermal  release  of  the  trapped 
charge.  Let  l,a  represent  the  quantum  numbers  of  the  trapped  electron 
and  positive  atom,  let  be  the  quantum  numbers  for  the  electron  and 

atom  after  the  electron  has  escaped.  Then  the  transition  or  release  rate 
u(la  -*  JP)  as  given  by  oi’dinary  time  dependent  perturbation  theory  is 

«(l,cc  -  J,3)  -  f  1  p  (A-T) 

where  is  the  density  of  final  states  and  where  the  energy  of  the 

final  state  must  equal  that  of  the  initial  state. 

The  actual  transition  rate  for  the  release  of  trapped  electrons  Is 
the  sum  of  w(iQ!  -  jp)  over  all  initial  states  weighted  by  the  Boltzman 
factor  and  over  all  final  states  consistent  with  the  requirement  =  E 

Kubo  has  calculated  this  sum.  To  do  so,  he  assumes  that  the  motion  of 
the  positive  atom  is  that  of  a  three-dimensional  harmonic  oscillator 
with  frequency  w  if  an  electron  is  trapped,  u'  if  not.  The  initial 
electron  state  is  hydrogenlike,  the  final  state  is  represented  by  a  plane 
wave.  With  these  approximations,  he  finds  that  the  actual  transition 


rate  Is 


W  = 


exp 


e 

lsl(a8  -  8®). 


^  large  (A-8) 

~  small  (A-9) 


In  these  formulae,  e  is  the  energy  difference  between  the  trapped  and 
free  states  of  the  electron,  8  is  defined  by  the  relation  (u'/u)  =  1 
A  =  256(m/M)((oa)' ,  where  m  and  M  are  the  masses  of  the  electron 
and  atom,  respectively. 


95 


APPMDIX  IX 


The  first  program  which  appears  below  was  employed  in  computing 
glow  curves  -  Eq.  (2-26)  -  for  L-phenylalanine .  These  curves  are  shown 
in  Fig.  20. 

The  second  program  was  used  in  computing  the  glow  curve,  Eq.  (2-65), 
for  a  four- unit  polymer. 


$$$  TIME  ON  1116 
WEDNESDAY,  6/6/62 

S538  JOB  6MIN  gill,  JIM  BIOPHYSICS  CONTACT  TESLER 

BURROUGHS  ALGEBRAIC  COMPILER  -  STANFORD  VERSION  7/25/61$ 
INTEGER  A,K$ 

ARRAY  11(5),  THErA(5)  =  (l.0T**3,  1.27**3, 

1.4t**3»  1.27**3,  1.27«*3),  gamma  (5)  =(l.25«*4, 

T.57**4,  k.h6**3,  1.25**4,4.46**5)$ 

WRITE  ($$TIT,LES)$ 

WRITE  ($$ISS,UES)$ 

LB=86$  RI=14 

FOR  T  =  (0.02,  0.02,  0.54)$ 

ENTER  INTEG$ 

LE=T3.06$  RI=2T.97$ 

FOR  T  =  (0.56,  0.02,  2.0)$ 

ENTER  INTEG$ 

LE=76$  RI-26.5$ 

FOR  T  =  (2.02,  0.02,  3.0)$ 

ENTER  INTEG$ 
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lb-8'3. BI-214 

FOR  T  -  (3.0a,  0.02,  4.0)$ 

BNTER 

WRITB($iG0Y)$ 

00  FLU3HIT$ 

SUBROUTINE  IHTBQ$ 

BEGIN 

te=i/(le+ri.t)$ 

TEA=!l/  (LB+HI  .  (T-0  .005)  )$ 

TBB=l/ (LE+RI . (T-0 .01) )$ 

TEC=1/(LE+RI. (T-0.015))$ 

TED=l/ (LE+RI . (T-0 .02) )$ 

FOR  A  =  (1,1,5)$ 

BEGIN 

THAT=-THETA(A)$ 

TEETH=TEBTH+ (TEETH-EXP (TEC . THAT ) +EXP (TEA . THAT ) ) $ 

bath=exp(teb.that)$ 

II  (A)  =  II  ( A) +(1 . 6666666**-3 )  (TEBTH-tTEEMH+BATH+BATH+EXP  (TED  .THAT ) 

+EXP(TE .THAT) )$ 

END$ 

WRITE  ($$VAL,IIES)$ 

RETURN 

END$ 

OUTPUT  VAL(T,  FOR  K=(l,l,5)$  GAMMA(k) .EXP(-THETA(K) .TE) .EXP(-GAMMA(k) 

.II(K))), 

ISS(T,  FOR  K-(l,l,5)$  GAMMA(K).EXP(-THETA(K)/86.0)), 

TIT(FOR  K  «  (1,1,5)$  GAMMA(K),  FOR  K=(l,l,5)$  THBTA(K))$ 


9T 


i- 


mms!  UB8(X4.2>  3(B2,  X18.8),  VO), 

I^(B2,  ^tQAIWA  »  *,  Bl,  3(B6,  F9.3,  Bj),  WO,  B2,  mSSSA  •  *, 
Bl,  3(b6,F9.3,B5),  W2,  Bl,  •»!-*,  WO,  Bl,  *—*,  W2), 
C0Y(B2,  #N0  MOBE  VALUES  OP  I(t)  TO  COMPUTE*,  WO)$ 

FLUSHIT..  FmSH$ 

COMPILED  PROGRAM  ENDS  AT  0552 
PROGRAM  VARIABLES  BEGIN  AT  9281 


TIME  ON  1109 
WEDNESDAY,  6/6/62 

S538  JOB  'fMIN  GILL,JIM  BIOPHYSICS  CONTACT  TBSLBR 

BURROUGHS  ALGEBRAIC  COMPILER  -  STANFORD  VERSION  7/25/61$ 
GAMMA=T.57**^ 

THETA-1. 27**3$  II«0$  THAT«-le7**3$SQ=SQRT(P)$  STaSQRT(3)$ 
SQA-(1+SQ)/1).$ 

S(?B=(1-SQ)/4$ 

sq,c»(i/sq-i)gamma$ 

S(5D=  ( -i/sq-i)gamma$ 

STA-(34ST)/12$ 

STB=(3-ST)/12$ 

STC-(ST/2-1)GAMMA$ 

STD- ( -ST/2-1  )GAMM4 
WRITE($$TIT,LES)$ 

WRITE  ($$I3S,  UES)$ 
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le»86$  ri-^ 

FOR  T  -  (0.02,  0.02,  0.5^*)$ 

ENTER  DfCEG^ 

LE=T3.06$  RI-27.97^ 

FOR  T  =  (0.56,  0.02,  2.4$ 

ERTER  Iin?ECr$ 

le=76$  ri  =  26.5$ 

FOR  T  =  (2.02,  0.02,  3.0)$ 

EirrER  IRTEG$ 

LE=83.5$  RI=24$ 

FOR  T  =  (3.02,  0.02,  4.0)$ 

EHTER  IRTEG$ 

WRITE  ($$COY)$ 

GO  FLUSHIT$ 

SUBHOOTINE  INTEG^ 

BEGIN 

TEE?rH=TEBTH+(TEETH=EXP(THAT/(LE+RI.  (T-0.015)))+EXP(THAr/ (LE+RI 
.(T-0.005))))$ 

BA!rHBEXP(THAT/(LE+RI,  (T-O .01) )  )$ 

II=II+(l  .6666666**-3)  (TEErTH-tTEETH+BATH+BATH+EXP  (THAT/ (LE+RI 
. (T-0 .02) ) ) +EXP(THAT/ (LB+RI .T) ) )$ 

WRITE ($$VAL, UBS )$ 

RETURN 

ENB$ 

COMMENT  THE  6  IN  THE  DENOMINATORS  IS  1/N(0)  $ 


<*•* 


OUTPUT  Via/(T,  (aAMMA*BXP(THAT/(US-t«I.T)))(BXP( -GAMMA. II) 

.EXP(SQC . II) .EXP (SQD . II )  +6TA .  BXP(STC . Il) 
4STB.EXP(STD.II))/6), 

I3S(T,  (aAMMA.EXP(THAT/86))  (l+SQA-t3qB4QTA+STB)/6), 
TIT(QAMMA,  THETA)$ 
format  UE3(x4.2,B7,X18.8w), 

LES  (B2,^^AMMA  =  *,BT,F9.3,W,B2,^‘THETA  =  *,BT,F9.3,W2, 
B1,^(T=*,W,B1,*~*,W2), 

coy(b2,  *no  more  values  of  i(t)  to  compute*,  W0)$ 

FLUSHIT..  FINISH^ 

COMPILED  PROGRAM  ENDS  AT  O56I 
PROGRAM  VARIABLES  BEGIN  AT  9225 
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